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I. SUMMARY 

The analysis of the acoustic liner which was presented in the 
Sixth Yearly Report has been extended to higher order and the effect of 
chamber flow is considered. The analysis is divided into two parts: 
one for a near-resonant situation and another away- from resonance. The 
second-order analysis of the latter part is presented in detail in this 
report (Section III-A) . Also a short discussion of the approach used in 
the analysis of the near-resonant case is presented. The chamber flow, 
including both the mean flow and oscillatory flow is shown to modify the 
equations which determine the cavity pressure and orifice velocity oscil- 
lations. Calculations are presently being performed so that the precise 
effect of the chamber flow as well as the cavity orifice geometry will 
be determined. 

Based on a numerical approach, the acoustic admittance has been 
calculated for a number of liner design conditions in an attempt to find 
possible simplifications in the representation of the flow field by this 
means (Section III-B) . The merits of this approach, as well as the prob- 
lems involved, are discussed. 

A series of tests performed in a rocket chamber environment to 
determine the optimum circumferential placement of acoustic dampers is 
described in Section III-C. Because this hot environment proved unaccept- 
able as a test bed, a cold flow approach is now being pursued and is 
proving quite fruitful. 

In another cold flow simulation study of liner behavior, shock- 
type waves are being generated and damped (Section III-D) . This study is 
providing insight to not only the damping efficiency with regard to non- 
linear waves; but, via use of shadow photography, is providing a visual 
history of the events that take place in the cavity e.g., the jet-like 
flow, vortices and turbulent phenomena. 

In the energy source study of steady and unsteady combustion 
based upon the use of a shock wave passing through the combustion zone, 
two sets of equations have been generated for the steady-state (Section 
IV-A) . Although complete experimental data have not been collected, 



preliminary experimental data falls in the expected calculated range. 
Further experimental results are currently being obtained to allow 
a complete solution of the steady-state problem. The unsteady problem 
now is receiving increased attention. 

In another area dealing with energy sources, the effect of 
tangential spin direction is being investigated as a chemical kinetic 
effect (Section IV-B) . LOX/liquid methane testing is described which 
yields interesting extensions to the previous LOX/ALC and LOX/RP-1 data. 
Reaction rate data from associated experiments (gas rocket and kinetic 
flow reactor) are also discussed in relation to this problem. 

Several analyses of nonlinear transverse oscillations with 
shock waves in annular chambers have been performed. The Priem-Heidmann 
vaporization model provides the combustion rate-controlling mechanism 
in the theoretical model. Droplet drag is neglected. An important ratio 
is that of the transverse velocity perturbation to the relative steady- 
state axial velocity between droplet and gas. In a paper to be presented 
at the Twelfth International Combustion Symposium by Professor Crocco, 
it has been shown that for small absolute values of this ratio, insta- 
bility can only be possible when triggered by finite amplitude distur- 
bances. However, no time-lag effect occurs due to the particular order- 
of-magnitude analysis which has been employed. In Section V-B, progress 
on another perturbation scheme which retains the time-lag effect to the 
order of interest is presented. 

Also considered in Section V-A is the case when the transverse 
velocity perturbation is large compared to the steady-state relative 
velocity. Here, it is shown that all oscillations decay with time. 
However, the analysis becomes invalid when the velocity ratio is of order 
unity. The small velocity ratio and large velocity ratio results combine 
to indicate that a stable periodic solution can exist when the velocity 
ratio is of order unity. 
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II. INTRODUCTION 

This is the Eighth Yearly Report on the research at Princeton 
University on the problem of combustion instability in liquid propellant 
rocket motors. These investigations are being sponsored by the National 
Aeronautics and Space Administration under NASA Contract NASr-217 and, in 
the past, under NASA Grant NsG- 99-60. The nonlinear aspects of the problem 
have been emphasized in both the theoretical and experimental portions of 
the program. Presently, the theoretical research is concentrating on 
studies of possible mechanisms of combustion instability such as droplet 
burning, the nonlinear gasdynamic oscillations associated with combustion 
instability involving shock- type waves in both the tangential and longi- 
tudinal modes, the effect of damping devices on stability with particular 
emphasis on nonlinear effects associated with acoustic liners, and funda- 
mental studies of the basic energy sources. In connection with the theo- 
retical program, an experimental program of equally broad range is being 
conducted. This program includes basic experiments with rocket motors and 
"cold" chambers which have been designed to study various phases of the 
combustion process and/or the instability phenomena to provide understanding 
which would prove useful in the construction of theoretical models. It also 
includes studies of the effects of exhaust nozzles and damping devices 
upon the instability as well as controlled testing of liquid rocket motors 
with parametric variations intended to determine the factors which strongly 
influence the incidence of nonlinear combustion instability. 

In the theoretical study of mechanisms of combustion instability, 
the greatest difficulty involves the construction of the model or models 
which contain the essential physical aspects of the phenomena. In this 
respect the basic experiments are extremely useful. The theoretical study 
of the nonlinear gasdynamic oscillations presents another type of difficulty. 
Generally, mathematical techniques for the treatment of nonlinear partial 
differential equations are not well- developed and, in most cases, they are 
undeveloped. Therefore, the theoretician must develop techniques as he 
works with the physical problem. 



6 


As is the norm for any status report, a considerable amount of 
the material discussed herein is in the state of active investigation. 

Thus, only tentative conclusions or observations are sometimes possible. 

At least a mention of one example of each type of work being conducted 
has been made so that this report would prove to be of maximum benefit 
to those readers working in similar areas. 

During the past year the following technical reports, papers and 
book have been published: 

"Axial Mode Shockwave Combustion Instability in Liquid Propellant 
Rocket Engines", C. E. Mitchell, NASA CR- 72259 (Princeton 
University AMS TR 798), July 1967. 

"A Preliminary Study of the Effects of Vaporization and Transverse 
Oscillations of Jet Breakup", J, A. Newman, NASA CR- 72258 
(Princeton University AMS TR 785), July 1967. 

"Nonlinear Wave Phenomena in Liquid Rocket Combustion Chambers", 

L, Crocco, C. E. Mitchell and W. A. Sirignano, presented at the 
Colloquia on Gas Dynamics of Explosions, Brussels, Belgium, 
September 1967. 

"Behavior of Supercritical Nozzles Under Three-Dimensional 
Oscillatory Conditions", L. Crocco and W. A. Sirignano, AGARDograph 
117, 1967. 

"Flow Behavior With Acoustic Liners", D. T. Harrje, W. A. Stinger 
and W. A. Sirignano, ICRPG 4th Combustion Conference, CPIA 
Publication No. 162, December 1967, pg 103. 

"Nonlinear Transversal Instability with a Simplified Droplet 
Evaporization Model", L, Crocco, Ibid, pg. 159. 

"Axial and Transversal Instability Using the n,^C Model", 

C. E. Mitchell, L. Crocco, Ibid, pg. 169. 

"Simulation of Droplet Vaporization and Jet Breakup in a Tangential 
Mode Environment", J. A. Newman, G. L. Waring and D. T. Harrje, 
Ibid, pg. 249. 

"Mixture Ratio Survey Technique Based on Propellant Temperature 
Differences", G. R. Black and D. T. Harrje, Ibid, pg, 281, 

"Experimental Investigations of Liquid Propellant Combustion 
Processes Using Streak Photography", R. M, Williams, Jr. , 

NASA CR-72371 (Princeton University AMS TR 800), February 1968. 
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In addition, two papers relating to this research will be 
presented at the 12th International Symposium on Combustion at Poitiers, 
France, July 15-19, 1968. 

"Research in Combustion Instability in Liquid Propellant Rocket 

Motors", L. Crocco. 

"A Theory of Axial-Mode Shock Wave Oscillations in a Solid 

Rocket Combustor", W. A. Sirignano 

Additional background and history of the research on combustion 
instability at Princeton University may be found in the partial list of 
References 1-12. In addition. Reference 13 represents the final report 
in an extensive series of reports on earlier work primarily concerned with 
the linear type of combustion instability (i.e., spontaneous buildup from 
steady- state operation). 

It should be acknowledged that this work made use of computer 
facilities supported in part by National Science Foundation Grant NSF-GP579. 
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III. DAMPING DEVICES 

A. THEORETICAL STUDIES OF ACOUSTIC LINERS 

In this section, an analytical study of the flow field asso- 
ciated with an acoustic liner is presented. This particular theoret- 
ical approach had been initiated by Sirignano, whose work can be found 

2 

in this group* s Sixth Yearly Report, The model used in the study 
presented here is much the same as the one developed by Sirignano. 

The basic assumptions used in the previous work are made use of here 
except that the assumption of no flow in the combustion chamber is now 
relaxed, viscous effects are omitted, and the effects of all harmonics 
in the chamber pressure are now considered. The assumed representation 
of the chamber conditions is done in such a way so as to represent any 
combustion instability type commonly experienced in liquid propellant 
rocket thrust chambers. In considering a chamber flow within this model, 
a clearer insight can be gained as to its effects on dissipation than 
that gained by the current empirical approach. A computer program has 
been developed to perform calculations in an attempt to correlate the 
theoretical predictions with experimental studies; however, the results 
of these calculations are inconclusive at this time and will not be 
presented here. 

Model 

Aside from the detailed justifications which could be found 
in the work by Sirignano, the basic description of the assumed flow field 
is as follows. 

Each orifice (hole) in the liner is considered separately along 
with its own individual backing volume. The system then geometrically 
constitutes a Helmholtz resonator. The low amplitude regime of potential 
flow is neglected so that there will exist a jet on the side(s) of the 
orifice in which flow is away from the orifice. No mean flow exists in 
the orifice, but both a mean flow and perturbed flow exist within the 
combustion chamber. For flow into the orifice, the flow just outside 
the orifice is considered to be quasi-steady. On the chamber side, when 



the flow is away from the orifice, the pressure that the jet experiences 
will be some average pressure (defined in terms of a pressure coefficient) 
which is not necessarily the chamber static or stagnation values because 
of the presence of flow within the chamber passing over this jet. When 
the flow is into the resonator cavity (backing volume) , it is assumed that 
the pressure across the jet at exit is equal to the static pressure within 
the cavity. Within the orifice itself, it is assumed that the flow is 
one-dimensional, isentropic, and unchoked. 

With this assumed model, the mechanisms for dissipation will 
be the following: 

1) the jet breakup produced on both sides of the orifice 

2) the eddies produced in the chamber flow due to the 
interaction of the chamber flow with the jet. 

The kinetic energy lost in the jet can be calculated directly from the 
analytical results. The interaction between the jet on the chamber side 
and the chamber flow produces a more difficult problem and a solution has 
not yet been attempted in connection with this work. It would be expected, 
however, that the second mechanism is not as important as the first since 
it occurs only during half of the cycle and since the jet velocities 
usually are higher than the chamber velocities. 

Formulation 

The mathematical formulation of the above model is now presented. 
The nomenclature used is defined in the table found at the end of this 
section. All quantities are non-dimensional ized with respect to the mean 
chamber state properties, the mean chamber speed of sound, and the orifice 
length. The problem will be first cast into a form convenient for the 
calculation of orifice velocity and orifice density. 

The governing equations for the orifice flow are the continuity 
and momentum equations with the isentropic relation replacing the energy 
equation. Thus, the following two equations describe the flow within the 
orifice 


P t + = o 


(l) 



For flow into the orifice on the chamber side, 


where 



n 






\rj 

*■ c P 

uf- 

z c p 7T 


= />' v‘/p^]/[ti-(^)v//^J 


Now from the isentropic relation 

P = P, ' 


so that 


the above could be written 



wr_ i*-' 

w x ‘ J 


For the cavity side, when there is flow into the orifice 


where 


4 = 


= Tr*= r' = ^^74 


>2 


f / - wp ~) ^ 


*F 




Since 


e, = >2 


4 


there results 
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For the cavity side when there is flow away from the orifice 

Fk ^ & 

For flow away from the orifice on the chamber side, some 
average exit pressure must be used as mentioned before. This is done 
by making use of a pressure coefficient, defined as (in dimensional form) 

0 >,» = (P, 

where P 0 is the static pressure at a point on the surface of the jet. 
Since the surface of the jet will be ill-defined and the jet will be 
distorted, this pressure coefficient will not in general be the same as 
the familiar one obtained for a circular solid cylinder in cross flow. 

It is also true that this coefficient will vary with time due to the 
unsteady nature of the flow field; however, since its effect will be 
seen to be of second order, for a second order result; both a time 
averaged and circumferentially averaged coefficient can be used to obtain 
the average pressure. Thus 

P = -tn. + Z 

Then for the chamber side with flow away from the orifice 

fi = LifpJi'/c, + gj* 

Equations (1) and (2) contain two unknowns whereas the above 
conditions in the cavity side contain an additional unknown, g. . The 
continuity equation written for the cavity is used to eliminate this 
additional variable. For the cavity then 




dLL Qtl 

dt 


( 6 ) 
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For flow from the orifice to the cavity, since 

6 =£*= P* 


the above becomes 


_ AJLo (A 

dt V r * 


(7) 


For flow from the cavity into the orifice, the isentropic 
relation (Eq. (4)) together with Eq* (6) yields 



( 8 ) 


The formulation of this problem then consists of the two 
governing relations (1) and (2) together with the above boundary condi- 
tions which are summarized below 



+ flow 


- flow 
+ flow 


- flow 


The boundary conditions at x = 0 illustrate the effect of 
the presence of a chamber flow on the orifice velocity. This effect is 
one of a modification of chamber static pressure, but this is accomplished 
in a different manner for the different flow directions. For flow into 
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the orifice, this effect is clearly a nonlinear one arising because 
of the modification of the flow field in general. For flow into the 
chamber, the effect is similar to an aerodynamic one modifying the 
local static pressure at the jet exit. 

Since the governing equations constitute a hyperbolic system, 
for the formulation to be well-posed mathematically, the boundary condi- 
tions must be open and specify both a value and a derivative at every 
point of the boundary (Cauchy conditions) . The boundary conditions are 
open in the sense of the time variable and the Cauchy conditions are like- 
wise achieved. This latter fact can be easily visualized by replacing 
LL by p through the use of the two governing equations and by noting 
that the condition at x = 0 supplies a value for p, (the chamber con- 
ditions assumed known) and the condition at x = 1 supplies a value for 
the derivative of ^ . The derivative of and the value for 
could then be found from the relation for p which must exist through 
the orifice. 

The method of solution will be one of small perturbations 
pursuing the periodic solution. In seeking a solution, a distinction 
must be made between the near-resonance and off-resonance conditions. 

Both solutions have been obtained, and the off -resonance case is now pre- 
sented in detail. 

Off-Resonance Solution 

The chamber static pressure and velocity are assumed to be 
given in the following form 

/? “ 1 + € * + €*£+••• 

l/ x = v" + eaUs* + e*a*V? } + 

where all quantities are assumed known with the proper ordering as to € , 
except for ; the mean flow, constant in time, which is assumed to be 

of perturbation magnitude. The quantities ^and Vx * > ( > 0 will 

be functions of time. The angle between the mean flow and the perturbed 
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flow is assumed known and is defined as 



For the case of longitudinal instability c< = 0° , and for the case of 
transverse instability o( = 90° . From the boundary condition, note 
that the basic solution (solution for 6 — ^ 0) will depend on the manner 
in which the limiting form is achieved; that is, it depends upon the 
direction of the flow when € — +* 0 . This, of course, is an inconsis- 
tency which occurs due to the simplifying assumptions made as to the 
discontinuous nature of the jet. A jet is defined in terms of the 
direction, with no regard as to the magnitude of the orifice velocity. 

It will be explained shortly how this inconsistency is treated. 

Substituting the above expansions into the W* functions 
yields formally: 

iVj = 2f + € F + £*S +OrC 6 3 ? 

where H , F , and S are various functions of Z//°* 3 P? * > and other 
similar terms. The above mentioned inconsistency can be treated so as 
to produce no mathematical difficulty by making use of the fact that l/^ o) 
is of perturbation magnitude and ordering with regard to this magnitude. 
When this is done here and in similar occurrences, the basic solution 
will formally reduce to the case in which the chamber static conditions 
prevail throughout the chamber, orifice, and cavity - for either flow 
direction. Thus, after neglecting small terms 

Wf = tffa y Z/f" . -L(*-l)al/?‘v?’cc:« , XjLSV?*] * OCe 3 ) 


so that, in effect, contributes only a second order term. One also 


obtains 
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The boundary conditions could now be expanded, and after 
neglecting small terms and equating like powers of £ 5 one obtains 
at x = 0 


for + flow 


e,"' = i 


e/" = 


- -L p 


O) 


- 4 & 


_ jfcZ 
^ y* 




+ 

<*) yr O) 


V"* - f < 




, I 7 r' a/ 7r l " 

+ -r <zUx cos 


+ — — 7 / 

* 6* Z V * 


(»r 


(9) 

( 10 ) 


(ID 


for - flow 



( 12 ) 



J~ p 
y <r 


O) 


(13) 




y -/ 
J1P 



_/ 

* 



/ 3? — yr^r r('> „ 
jrCpG-Vz ££ COS X 


/ 

6‘ 


J. n 1/ 

2 ^P ‘Or 




(14) 


at x = 1 

for + flow 



(15) 


( 16 ) 
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for - flow 



o 

( 17 ) 

< ir‘ = ° 

( 18 ) 


Note the occurrence of the term •— f or X s 2 ) 

in the second order boundary conditions* This term embodies the second 
order "jet effect" captured by Sirignano. It is interesting to note 
how this term corresponds to a jet effect. The origin of the term 
itself is due to the nonlinear effect of inflow to the orifice, but its 
significance is enhanced because of the absence of the term upon outflow 
(jet). If the outflow were also quasi-steady, this term would appear 
again. 

For the purposes of dissipation, the liner system is treated 
as a jet diffuser with a zero efficiency, that is, none of the kinetic 


energy of the jet is used in regaining stagnation pressure but is merely 
dissipated through viscous action. From this consideration, it is 
reasonable to postulate that the further the system is away from isentropic 
deceleration of the jet the better. Thus, when C p = 1, the exit pressure 
of the jet is the stagnation pressure in the chamber and this indicates 
a poor condition. This is illustrated by the second order boundary con- 
ditions at x = 0 , for when C P = 1, the equations become the same 
except for the term, and in this case, the action of the discon- 

tinuities in the equations for different flow directions as a 
forcing function is diminished. This interpretation is also conformable 
to the fact that the smaller the algebraic value of C p , the greater is 
the pressure difference across the orifice and thus the acceleration of 
the fluid. From the above reasoning, in practice it may be advantageous 
to shape the orifice cross sectional area so that the shape of the dis- 
torted jet enhances this aerodynamic effect. 

The relation between ^ and yg. is desired, and upon 


expansion: 
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for + flow 



p 0) _ J_ p 0) 

K Z ~ ^ rk 

± [ P ( *> _ u IzJ P ou 7 

Y l & 2 y <* J 


(19) 

( 20 ) 


for - flow 


p ( ° — jL p (,) 

\z ~ Y a 


( 21 ) 


c = tfe"- i ¥ z o> ‘] - i «r w2> 


The expansions of the governing equations become 


continuity 



+ U? 


= o 


(23) 


momentum 


+ r p (/ v /J ; x = o 
< /; + p* ,; = * 


(24) 


(25) 


* p"' * - (*-r)? 01 ?x - * 

First order solution 

Equations (23) and (25) yield the wave equation upon taking 
partial derivatives and equating mixed derivatives: 


U 


o) 


xx 


U 


o) 

it 


(27) 
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a*' 7 

xx 


Q<!) 

If* 


The most general harmonic solution periodic in time is of 


the form 


T Cx) e 


imt 


where J(x) may be complex (solution at different stations may differ 
in phase), and m is related to the frequency of oscillation. Then, 
assuming 

u m =[?,(*) / 

f‘" = + C^cx)] e‘ mt 

Equations (27) and (28) yield 



— 


+ Be"”* 



C e <M * 

+ oe e/r> * 


-r 

Se"”* 

+ Fe e '~* 


— 

Ge**”* 

f- Ke c ~* 


Either Eq. (23) or (25) could be used to determine the relations between 
the arbitrary constants. Doing so yields 

a a> * Cc.*0*** ) + cc Cc‘ m * + Oe "” x ?] e 6 ”* 

[(-Ae"”* + ) + c'C -Ce‘ M * +o<s- e/ ”*)3 e"** 

Expanding and collecting real parts gives 

5 sin m x] cos mt + [-Nc03/»X s/rnnt 


(28) 


U°' = [m cos mx 


(29) 
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f'" = [■* cosmx + Aisinmxjcosntt + [S cos m x + M sin sin/nt 

Note that the first order boundary conditions are the same 
for both + and - flow. Assuming — cosot and using the boundary 
condition at x = 0 to evaluate Equation (30) gives 


m — co j ft 


-L , 5=0 


(30) 


(31) 


The boundary condition at x = 1 together with Equations 
(29) and (30) yield 

[Rcoso ~ N s/o oJcj sin o t + cj M sin cj cos cat 
- (g[/icos co cos oit -( ft cos to + ftsin oj)sin<otJ = O 


Thus 


0 


.. _ f cocos co /• gsincj l a 
N / OSinui — ^eoj <*> J ° 


(32) 


or 


U {,) — yr[dcosa»x + sin ci x] sin at 
f‘" - il-A sz/fojx 4* cos <jj cos cjt 


(33) 


where 


A — C oc0$<*> +Qstn co 

^ Gj&frto* — pcosu* 


Now cj = £. nfj/c m 2 7rJP/A. 
so that for A »J CJ is small, and thus 

u M - tjp + 

f‘"“ i[- *j]cosat 
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Since the Helmholtz resonating frequency can be expressed as (for 
negligible end correction) 

or 

(3 = ( 2 7T f 0 J/c} 2 ' — Qjf 

for chamber oscillatory frequencies near this frequency, the above 
reduces further to (assuming approximate accuracy) 

u ‘" = y sina3t 

f"' = i[- + 1 


The last result is the first order result obtained by Sirignano, and 
note that in this case, the first order velocity is constant throughout 
the orifice with a linear density (or pressure) gradient 

Second Order Solution 


Evaluating Equations (24) and (26), taking partial derivatives. 


and equating mixed derivatives gives the inhomogeneous wave equation 
for *«V p“: 


~~ — [ M s/n 2 cox + A/cos£a)xJ s/ft2a>t 


_ = jj~ A/ s/n2cox + Af cos2coxjcos2a>t 

-f- P sift 2 ojx 4 Q cos 2 csx 


(34) 


tl - + yX^-1)] 

N - - -*jr[fO + r)] 
P= -&r[/>C3-r)] 

Q = fg[( 3 -r)(/-A‘>] 


where 
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These results are of the form 

= f, C£ CO x) s //)£&)? 

~ Px** — (2.00x7 + %(£ u>x7 cos 2 cot 

A particular solution periodic in time can be obtained by 
separation of variables and is 

l/p — [ B, s/r)£u>x + cos 2 cox - /V x s/s?2cox 

+ M xcos£oox]sJ»2cot 

^ C p } ~ 4 CO* [ Pc/r>2coX + QcosZcoxJ + B 3 X ~t B+ 

+ £B S Sfn2cox + BfCos 2cox — 

- AJ x c os 2 co xj cos 2 co t 

where B^ , - - - , B^ are to be determined such that this is a particular 
solution of Equations (24) and (26) . Substituting the above into either 
of Equations (24) and (26) will give the required result and is 

U 1 ^ = cos 2 cox “* Ns/n£ coxj x <s/n2 co t 

Ps/h£ cox 4- Qcos£cox] 

+ 4 ^J~[Cz<o P — M s/r>2cox + C jiz Q — A/x)cos2<ox]cos2-<*t 

where any arbitrary constants still remaining from the higher ordered 
equations have been set equal to zero. The complete solution will consist 
of the homogeneous solution and the above particular solution with the 
homogeneous solution adjusted so as to conform to both the particular 
solution and the boundary conditions. These boundary conditions will now 
be developed. 


■W'-' 
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The specification of the input functions will first be 

made 

oO <»* 

B CZ) ~ £ Cp n Cos*oit + A. %p jn 

i?~Z * 

IT X C> = e^coscot -h q u sintA>t 

The boundary conditions change for the different flow direc- 
tions. If the exact points of flow reversal were known (which in general 
would differ for each end of the orifice), then a Fourier expansion 
could be used to obtain a single expression for the boundary conditions 
valid for both flow directions. Since these points are not known until 
the solution has been obtained, one method of attack would be to treat 
these points as unknowns, solve the problem, and then evaluate these 
unknowns from the solution. A vast simplification can be obtained here 
by making use of the points in time when the first order velocity goes 
through zero. This is particularly convenient since the first order vel- 
ocity goes through zero uniformly in x . The justification of this 
approach leading to a result accurate to second order is the following. 

Let A be a point in time in which the first order velocity goes through 
zero. The actual velocity at t = A will then be of second order in 
magnitude. Since the slope of the velocity curve vs. time is of first 
order, the error incurred in the time of actual flow reversal will be of 
first order; that is, the actual time of reversal will be A 4* (7(€) , 

Now since this last quantity appears in the limits of integration necessary 
to evaluate the Fourier coefficients, the error incurred in the second 
order solution is of higher order so that the complete solution will be 
accurate to second order. 

Proceeding in the above manner, the boundary conditions at 
x = 0 are first written as 

- -j- <z*Vx** + 4" G ot 

+ e* -X Z '* + flow 

£Cuyf) = ic, + i-C.a t/J 0) t4'cos 

, J- r 2. r/ c * z 

€* t "P 2 Kr 



flow 
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The function to be expanded is then 




^ ~ (o»t) 




o < 77" 

7T< cot <Z7T j etc 


after performing the expansion, the boundary conditions at x — 0 
become: 


c = 


Go 4 G t COS cut 4 G p S/A<ot + G a cos /cot 

n ~ ^ cosAQjf , r r 2- £*±*££* 

4- Gj Stn /cot 4 G s ~~ZfT7P r £rz / ~ /7 ‘* 

/> n even 

+ GX + 6 e £ 7^7 

7 ft*/ & J?x/ *r ^ 

n *44 ff *44 

mm «►* 

+ G,Z- 7pT7 + <?,.£, TRWTft 

* e ve/* n 

* GXC „ cos n<ot + q^ s/s scot o 


(35) 


where the coefficients are defined as 

G 0 = i<?CeZ + 9 u)Cj 4- U P J 4 £ i VF'C* 4 

+ -fjf cl q u l^eos -?/>! - y4*cos £ co - 
G, ~ al/J^eu cos oc. ( i + £ p y 

G a = r T at/J\«cos~(j. 4 C F ) + y 2 r 

G 3 = w<z*CeZ + tfX*+e P ) 4 if*A*co**co - 

G4. ~ ^ G? G u <£ct C 1. + Cp ) 
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Gf — jj- q? $•* c i ~ Cp y 

G a = -£■ a trj*^ cos oc( J-C p ) 

G, - tti«*Cel+&) + £or? J *]CJ~Cp) 

&q — £^zr ^ G e u 2 % )C 1 ~ Cp ) 

G 9 — ~e 7r a Uz 0 ) e u cos « Cl - C p ) 

G /e> = ryZ A* cos* co J G„ — y 

As a check, it can be verified that for Cp = 1 , the 
boundary condition reduces to almost the same thing as that for flow 
into the orifice - the only difference being due to the y U C / )Z term. 

The boundary condition at x = 1 is treated in a similar 
fashion. The function to be expanded is 


O c*)t< 77 


X Ccot) = 


CO 

£)S‘ 


J^Acosgj Sin J] Sf/?JZu>'£ 77" c ot‘<^E7T > 


so that at x = 1 , after the expansion is performed 


C2) t 

=■ T'SwZoit + j* 2 C 4,-rp* cosnCjjt 

r? odJ 


( 36 ) 


where 


= ~j[ y> ^ [COS U> ~ds/»cJJ[/lc0S CO + S/rtto] ^ 

CaSOJ S/*Cj] 

Jz = yr^ [ ^ cos 00 * s/ " "7 
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The most general harmonic solution of the homogeneous version 
of Equations (24) and (26) is given by Equations (29) and (30) . The 
presence of constants and higher harmonics in the above boundary condi- 
tions indicates that for this homogeneous version, a solution with terms 
other than first harmonics is necessary. Thus, considering all possible 
solutions, it is known a priori that whatever the solution of the homo- 
geneous equations, it can be represented by 


l/**’ ~ COS POJX - Sjs/o/nox) cos/? cot 

+ C -/iGcos/rcox - fit, S//7/rcot] + E 

mo 

H [C ~ Eficossrtox NnSin/lctix) cos /teat 

H f7-t f 

+ C S^cos/zcox -f- s/yf/*eoy) S//r/fcot j + J) 


(37) 


Thus, at x = 0 , 

+ (f** = COS/? Utt ■+ S* sinneot] + 0 

+ 4 g?& + G) cos £ cot 

Matching the above with Equation (35) gives 

D » O 
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Making use of Equation (36) yields 

— ~ [ + C zzv G — s/fiZcot 
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All of the constants in the homogeneous solution have been 
evluated so that the second order solutions for U. ta> and Q Ci> are now 
complete. 


Cavity Pressure 

The cavity pressure could now be obtained by making use of 
Equations (19) to (22). In the second order expressions, a Fourier 
expansions must be done. The function to be expanded is 


%Cgj t) 



cos co + sin co 


C COsZcot — 1 ) 
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The results of this expansion give 
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Resonance Solution 

The details of the near-resonance solution will not be 
given here, merely the highlights of this solution will be presented. 

The characteristic difference of this solution from the above (from which 
all other differences follow) is that the orifice velocity and density 
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fluctuations are ordered so that they vary as the square root of the 
chamber pressure fluctuations. With this ordering then, the pressure 
work done on the system will be balanced by the kinetic energy dispersed 
by the jet. 

The algebraic manipulations become exceedingly more involved 
- the difficulty being the fact that an (/j + 1) order study must be 
made in order to obtain a complete (/?) order result. Because of this 
difficulty, it is not permissible to expand the boundary conditions 
based on the values in time in which the first order velocity goes 
through zero; namely because second order equations must be used in order 
to determine the first order velocity, and if the second order equations 
are inaccurate, this will reflect upon the first order terms. A correct 
method of attack was alluded to previously; i.e., the points of flow 
reversal are treated as unknowns in the problem, and once the final 
expressions have been found, these unknowns can be evaluated from their 
very definition. 

A more detailed explanation is as follows. It is assumed that 
at the ends of the orifice, the period of velocity oscillations is div- 
ided into two regions: 1) a jet exists on the cavity side (+ flow) and 

2) a jet exists on the chamber side (- flow) - see Fig. 1 . Let (p, 
be the value of a> t for which the velocity at x = 0 goes through zero 
and in which positive flow exists immediately after this point. Let 
9>/+ P, 8 be the next value of cot for which the velocity passes through 
zero. Negative velocity will then occur immediately after this second 
point, and the next point of flow reversal will be cot * + 2.7T - the 

8 

start of periodicity. At x = 1 , the corresponding values of Pa and 
are similarly defined (which will in general differ from corresponding 
values of p, and (p t ) . The solution of the problem is then formally 
with the (p ' s , and an expression for the velocity accurate to second 
order can be obtained such as 
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The p's can then be found from the following four algebraic equations: 

u(o, <P* 3 0 * 0 * <?/) = o 

c< (o, pa p e , PA AA ) 

«Ctj = o 

u cj- , ?/, AA &*) 

The solution in which positive velocity occurs after and 9 is 

desired. Note that any solution of the form = fy n ^ 7r is accept- 
able since the <p '*S appear only in the limits of the integrals necessary 
for evaluating Fourier coefficients. The justification for using the 
value of cof in which the approximate velocity (accuracy of (3^(e*)) 
passes through zero follows the same reasoning justifying the procedure 
used in the off-resonance solution. 

It should be noted that the above four equations cannot be 
solved as a unit. Four more equations exist for the determination of the 
arbitrary constants of the various homogeneous solutions - two constants 
for first order (from second order equations), and two constants for 
second order (from third order equations). Thus, for calculation purposes, 
it will be necessary to numerically solve a system of 8 algebraic trans- 
cendental equations . 

The region of validity of this near resonance solution is 
dictated by the expression 

[as tan CO - ^ J = CO) (38) 

In applications, $ is much smaller than unity so that this rules out o> 
becoming very large, and thus tano> will be greater than zero. It is 
interesting to note that from this criterion, the solution can be valid 
also for cases in which the chamber frequency of oscillations is far 
away (although smaller than) the Helmholtz frequency. It is also inter- 
esting to note that this is the criterion which predicts when the off- 
resonance solution becomes invalid (see Equation (33)), Thus, the two 
solutions are complementary. 


=- O 
= O 
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The details of this solution have been performed and a 
solution accurate to order e d has been obtained. 

Future Work 

Having obtained the analytical results, the point is reached 
where significant feedback can be made between theory and experiment. 

One direction of investigation will be to determine what geometrical 
configurations of the liner system lead to optimum dissipation. Certain 
trends should be predictable from these analytical results, and the 
facilities here at Princeton will allow both cold flow and hot flow 
experimental verification (the experimental programs are discussed in 
other sections of this report). 


NOMENCLATURE 


d proportionality constant relating chamber pressure 

fluctuations to chamber velocity fluctuations 

A orifice cross sectional area 


Cp 

Cr 

c P>n 


f 

2p,p 

2 * 

J 

p 

t 

Vz 

u 

X 


average speed of sound 

specific heat at constant pressure 

average pressure coefficient for jet on chamber side 

th 

coefficient of the n harmonic (cosine) in the 
second order chamber pressure 

coefficient of the cosine component in the first order 
chamber velocity 

frequency of chamber oscillations 

th 

coefficient of the n harmonic (sine) in the 
second order chamber pressure 

coefficient of the sine component in the first order 
chamber velocity 

orifice length 

pressure nondimens ionali zed with respect to the 
mean chamber static pressure 

nondimens ional time = (actual time) (c/jf) 

chamber velocity nondimensionalized with respect to c 

orifice velocity nondimensionalized with respect to c 

distance along orifice measured from chamber side divided by JL 
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NOMENCLATURE - continued 

]/ cavity volume 

oc angle between chamber mean flow and chamber perturbed flow 

Y ratio of specific heats 

$ nondimensional pressure amplitude parameter 

Co nondimensional frequency of chamber oscillations = 2 7rfJ/c 

Subscripts : 

& stagnation conditions in chamber 

1 conditions at chamber end of orifice 

2 conditions at cavity end of orifice 

I static conditions in chamber 

XI conditions in cavity 

Superscripts : 

0, 1,2, . , . indicates order of term 

positive flow direction is from chamber to cavity. 

nondimensional izat ion indicated above does not apply if noted. 

above subscripts do not apply to constants and incidental 
functions within solution details. 

B. LINER ACOUSTIC ADMITTANCE 

The admittance of an acoustic cavity is defined as the ratio 
of the velocity to the perturbation pressure in the chamber of the rocket. 
If both these quantities are sinusoidal, then the admittance will be a 
complex number, carrying magnitude and phase information. If the velocity 
in the orifice and the chamber pressure are non-sinusoidal, then there 
is a separate admittance for each Fourier component of these quantities. 

The investigation attempts to find out what the admittances are for a wide 
range of conditions, and use the information gathered to try to define a 
single, useful admittance coefficient. 

The computer program written by Stinger, and described in the 
last yearly report, was used to solve for the flow conditions. The assump- 
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tions made to arrive at this formulation were that the velocity in the 
orifice is constant at any time, and the orifice flow is isentropic, 
with different mean temperatures in the chamber and the cavity. A 
complete derivation of the program can be found in the last yearly report. 

The input variables are: The chamber pressure, in its 

Fourier components; $ = Ajf/V , which is an inverted measure of the cavity 
volume; ©* , which is a measure of the cavity temperature; A1 , which is 
the perturbation chamber pressure magnitude divided by the mean chamber 
pressure; and FREQ, which is the frequency of the oscillations divided 
by the Helmholtz natural frequency. 

Two different wave shapes were used as input for the chamber 
pressure, one of them sinusoidal and the other a typical wave found 
experimentally. For each of these wave shapes, the amplitude was varied 
from .1 to .4 times the mean chamber pressure, o< was varied from .05 to 
1.0 in six steps, and FREQ was varied from .6 to 2.2 . Q was kept con- 
stant at the value of .142 corresponding to the test rocket at Princeton. 
This was done mainly to keep the volume of data down, and partially 
because the effect of small changes did not seem to be significant. If 
further investigation into this problem is deemed desirable, one of the 
possible directions is to vary (3 

Four sets of graphs have been obtained in this study. The first 
two sets present the ratio of the second frequency velocity component 
divided by the first frequency velocity component for each of the two 
wave shapes. This velocity ratio will indicate when it is necessary to 
consider the second frequency component of the admittance ratio. The 
other two sets of graphs present the admittance ratio versus the frequency 
ratio for each of the two wave shapes. For the sinusoidal case, both the 
real imaginary components are plotted with the magnitude of the admittance 
ratio. For the actual wave shape, just the magnitude is plotted for each 
of the two frequency components of the admittance ratio. If the real and 
imaginary components are also desired, they are available. Typical 
examples of these various plots are presented in Figures 2 to 7 of 
this report. 

In general, the velocity ratio for the sinusoidal case has the 
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same overall shape for all the graphs, with an increase in the frequency, 
or a decrease in either Alpha or the amplitude causing an increase in 
this ratio. The velocity ratio also generally hovers around 0.2 for the 
Alpha's under 0.75, which is a normal running condition. Figs. 2 and 3 
illustrate these trends. This is a good indication of the nonlinearity 
of the problem, since for this case there is no second frequency component 
of chamber pressure, giving an infinite admittance. 

For the case of the actual pressure wave shape, the average 
value of the velocity ratio is generally higher, as would be expected with 
the inclusion of a driving pressure at the second frequency. Again, the 
shapes of the velocity ratio curves are quite similar over the range of 
Alpha and Al. However, this shape is quite different from the shape found 
in the perfect sinusoidal case, even though the second frequency component 
of the chamber pressure is only .18 times the first frequency component. 

The principal change in shape is a significant dip around the resonant 
frequency (see Fig. 4 ). Since at this point the velocity ratio goes 

to about ,1, it would be a good approximation to just use the first fre- 
quency component of the admittance around resonance. For large values of 
Alpha, the velocity ratio is almost constantly under .1, and the inclusion 
of only the first component of the admittance ratio could safely become a 
general rule. For these cases, however, Alpha is .75 or higher, which is 
abnormally hot for the cavity. 

The admittance for the first frequency component seems to be a 
fairly "nice” function for the perfect sinusoid case, especially for high 
amplitude. When the amplitude of the pressure is .4, the admittance is 
fairly constant. It should be noted, however, that for the cases with 
low values of Alpha, the flow became choked with the pressure amplitude 
of .4. For the case of the amplitude equal to .25, the admittance is also 
quite flat, centering around 2, for all values of Alpha tested (see Fig. 5 ). 
For the case of the amplitude equal to .1, the admittance varies from 
around 2 to 3.5. 

Although the first reaction to these results may be favorable, 
it can soon be seen that they are not so pleasant. The major change in 
admittance can be seen to come from changes in amplitude of the pressure 
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waves (compare Fig. 5 and Fig. 6 ). Since in the process of damping 
out a disturbance the amplitude changes continually, this sensitivity 
makes it more difficult to define a single admittance. 

A comparison of the actual wave admittances with the perfect 
sinusoidal described above show that the first component of the admit- 
tance for the actual wave behaves the same way that the perfect sinusoid 
does. However, the second component of the admittance has a different 
shape. The main difference is a dip around resonant frequency, that is 
a result of the decrease in the second component of velocity (see Fig. 7 ). 
As for the sinusoidal case, the curves are more constant for higher 
amplitude. For these cases, it appears that a reasonable approximation 
is to say that the second frequency component is zero around resonance 
(FREQ = 1.0 + .2), and equal to the first component of the admittance at 
other frequencies. For the low amplitude case, the second frequency com- 
ponent of the admittance is also more erratic. 

The accuracy desired will, of course, have to determine how 
far it is deemed reasonable to go with these approximations. There is 
some question as to the validity of the underlying assumptions of this 
analysis, so experimental verification will definitely be necessary. 

If the predictions of this analysis are borne out, then the admittance 
can be obtained from a graph, or a general value such as 2 can be used. 
The ability to neglect the second component around resonance and treat 
it as equal to the first component away from resonance is a good simpli- 
fication, but it ignores the main complication which is the great sensi- 
tivity of the admittance to changes in amplitude. This is especially 
critical when the chamber pressure will be going through a series of 
amplitudes. In the case where the amplitude is constant, the admittance 
can be fairly well defined. In the case where there are reasonable 
changes in the amplitude of the chamber pressure perturbation, then an 
average will have to be taken over a large range (from typical results 
as shown in Figs. 5, 6 and 7) or the idea of a single admittance will 
have to be discarded. 
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Figure 4 


Frequency ratio 






Admittance ratio 
a = . !95 Al = . 25 



Frequency ratio 




Admittance ratio 


- 40 - 



Figure 6 
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C. PLACEMENT OF DAMPING DEVICES 

Often in the cases of practical rocket motor design there are 
only limited locations where damping devices can be placed. More spec- 
ifically, for liner-type damping approaches, it would prove quite useful 
to know the most effective axial and circumferential locations for the 
placement of acoustic dampers. For the axial locations, tests at P and W, 
Lewis Lab - NASA, RPL and Princeton have all shown the location closest 
to the injector face to be far more effective than locations downstream. 
With regard to circumferential locations, far less is known. However, 

placement in limited circumferential locations, e . g. , at 120° intervals, 

1 2 

has been previously reported by Princeton. 5 The purpose of the 

current effort is to study the case of the liner used in conjunction with 
a minimum baffle. The baffle fixes the mode (e . g. , a three-b laded baffle 

would still allow a third tangential standing mode to exist between the 
blades) thus allowing the acoustic resonators to be tuned for a specific 
frequency. Since the effective damping range of these devices covers a 
very narrow frequency band, instability mode restrictions are highly 
desirable. Therefore, a prime consideration in this study was to deter- 
mine which circumferential locations for the acoustic dampers were the 
most effective. 

Tests were made on the 1 x 16 injector (sixteen spuds of like- 
on-like design), which experienced first tangential mode instability of 
250 psi pk-pk amplitude following a tangential pulse. These LOX/ ethanol 
tests were made to check the concept of using a limited number of acoustic 
dampers to bring dynamic stability to an otherwise nonlinearly unstable 
system (characterized at 150 psia chamber pressure). Damping was accom- 
plished with as few as three resonating cavities spaced at 120° intervals 
around the circumference of the liner. The acoustic damper had six, 

\ inch long, .100 inch diameter, orifices entering an acoustically tuned 

1,2 

cavity. It was found that the time to damp was short, being of the 

order of 25 milliseconds for nine active cavities and 17 milliseconds for 
27 active dampers. A cavity volume of approximately three times the 
design volume was reached before marginal damping was experienced. Design 
frequency of 607 o of that desired could still provide stable operation. 
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In summary, for this chamber-injector configuration, only a 
limited number of acoustic resonators were necessary and these provided 
damping over a sufficient band to eliminate the finite amplitude distur- 
bances encountered. 

Certain injector arrangements, where the propellant flow was 
made circumferentially non-uniform, caused standing modes which often 
avoided the locations where optimum damping was provided via the reso- 
nating cavities. Hence, baffles were used to control the modes and maxi- 
mize the damping. It was found that with such a three-bladed baffle the 
instability was damped in less than 10 cycles and hence, it was difficult 
to maintain instability for testing damper effectiveness. The mixture 
ratio variations (and hence temperature striations) were apparent from 
the marking on the aluminum oxide chamber wall coatings. Cavity measure- 
ments indicated an out of resonance condition immediately following the 
passage of the initial high amplitude wave. Actual conditions were some- 
times far from those assumed in the design. For example, in one case 
the cavity frequency was initially depressed more than 10% and required 
10 cycles to reach nominal frequency (2780 Hz). Beyond that point there 
was a phase lag of 138° across the orifice. If the phase map constructed 
from the nonlinear flow behavior model discussed in the Seventh Yearly 
Report'*' (NASA CR 72270, Figure 25) was used, the implications were that 
the cavity temperature was only of the order of 600°R, even after a regime 
phase reading had been reached. Another cavity in the same test had a 
phase lag of only 78° which would be representative of about 2000°R tem- 
perature. Tungsten-rhenium thermocouples placed around the perimeter 
showed there was a relatively cool region from the center of the spud, 
through the fuel portion and approximately 60% of the distance to the 
adjacent spud. In the cool region the temperature stabilized in the 500 
to 1000°R while in the hot zone it remained in the 2000-3000°R range. 

(These temperatures were in agreement with the phase measurements.) 

Along the zone interfaces there were severe temperature fluctuations. 

The spuds were reoriented in an attempt to equalize the wall 
surface temperatures by putting the fuel rich zone toward the wall. It 
was also hoped that the new configuration would lead to greater instability 
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because damping was occurring in about 50 milliseconds with only a few 
cavities active. Stability was still easily obtained thus presenting 
an unsuitable test condition. Therefore, an unlike-doublet injector 
was used since it was known to be inherently unstable. A limited-length 
cruciform baffle design to select the second tangential mode was used 
in conjunction with the liner. With the liner adjusted to the second 
tangential mode, a first tangential mode oscillation of low amplitude 
was observed in the chamber. When the cavities were tuned for the first 
tangential mode, the second tangential mode appeared in the chamber. 

A three-bladed baffle was used instead of the four-bladed design, and 
groups of cavities (using all the available liner) were tuned to eliminate 
the first and second tangential modes. This reduced the amplitude of 
both, but the first radial and several combined modes appeared. In 
summary, for this unlike injector the instability could not be damped 
with any combination of cavity tuning. 

A combination of spuds from the like-on-like, which was stable 
under all conditions, and from the unlike impingement, which was extremely 
unstable, was tried in an effort to reach a marginal instability range 
which was necessary to determine the optimum placement of the resonating 
cavities at pressure or velocity nodes. However, the combined injector 
scheme was still unstable and more closely resembled the unlike tests. 

Rather than continue in the quest for a proper rocket test bed 
for determining optimum circumferential locations for the acoustic dampers, 
attention was turned toward cold flow simulation. Perhaps one of the 
principal sources of difficulty in the hot firings was that the level of 
uniformity reached in higher thrust multi-element injector designs cannot 
be achieved in 12 or even 16 spud motors. This would appear to be partic- 
ularly true for the wall conditions where barrier coolant improves uni- 
formity for the full scale combustors. 

The cold flow testing is being conducted in an apparatus formerly 

14 

designed for oscillating heat transfer research. A new test section 

was fabricated which included a variable volume cavity and a connecting 
orifice equipped with hot-wire anemometer ports so that local flow can 
be measured. Pressure in the cavity and in the main duct (at the damper 
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location) are measured by flush-mounted Ristler 601A transducers. 

Hot-wire probe ports are also included in the duct. The test section 
is mounted on an axially traversing section of the duct so that both 
node and anti-node locations can be checked for duct longitudinal mode 
frequencies of 270 cps or higher. Because such longitudinal velocity 
and amplitude characteristics are very similar to the baffle cavity 
modes, this survey should provide information on optimum circumferential 
placement of the damping cavities. In addition, and possibly more im- 
portant, will be acquisition of pressure and velocity data with regard 
to both phase and amplitude. These data together with the nonlinear 
studies discussed next should provide an improved basis for optimum 
damper design. 

Preliminary results have indicated improved damper efficiencies 
with the volume in excess of the design volume. This is in agreement 
with the hot firing data previously mentioned, where lower frequency 
(higher volume) designs proved highly effective. The results of these 
studies will be reported at the 5th ICRPG Conference this fall. 

D. DAMPING OF SHOCK-TYPE WAVES 

One aspect of the problem of acoustic liner design is the 
effectiveness of damping shock-type waves. Since this type wave in the 
tangential mode is characteristic of some of the large diameter chambers 
as well as in some of the smaller diameter chambers with reactive propel- 
lants (JPL experience in the 11-inch chamber^ ~*) , and is representative 
of the longitudinal mode, this study was initiated to determine how a 
liner behaves under such an environment. 

The apparatus used in this experimental study (as shown in 
Fig. 8 ) consists of a shock wave generation device, a variable- length 

shock tube (with a windowed resonator cavity segment) , pressure trans- 
ducers and photographic equipment. The high amplitude shock wave is gen- 
erated at one end of the tube (referred to as the ''chamber 11 ) and rapidly 
steepens into a planar shock which subsequently interacts with a properly 
tuned Helmholtz resonator (referred to as the "cavity"). The following 
describes the preliminary, qualitative findings from a study of that 
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interaction. 

Data taken were of two forms, pressure vs. time plots and 
shadowgraph pictures. The pressure plots yielded information on damping 
rates and relationships between cavity /chamber wave amplitudes. The 
photographs portray a single cycle in detail, illustrating jets, vortices, 
and shock wave deformation. 

First, a more detailed description of the experimental appa- 
ratus is necessary. The chamber is 12 inches long and has a one inch 
square cross section. It is segmented, allowing the cavity segment 
(Fig. 9 ) to be placed almost anywhere along the chamber length. The 
cavity section has windows which allow observation of the entire cavity 
and that part of the chamber immediately adjacent to the orifice. The 
acoustical volume of the cavity can be continuously varied (via a square 
piston) and the orifice diameter can be changed from zero to three-eighths 
inches (via a removable orifice section) . 

The pulse gun consists of a spring-powered, cigar- shaped, ram 
which slides through the "hole 11 of a doughnut -shaped, high-pressure plenum 
chamber. The ram is hollow. A half- inch long port on the ram allows a 
burst of high pressure gas (N^ to enter the ram’s interior (and thus, the 
chamber) as the ram slides quickly by the plenum. The pulse generated by 
the pulse gun is actually two shock waves - a low amplitude wave immedi- 
ately followed by a much stronger wave. This pulse gun was used as a 
matter of convenience and was considered adequate for this preliminary 
s tudy . 

The pressure sensors used are Kistler 601 transducers. They 
are installed as shown in Fig. 9 . These signals are amplified, recorded 
on a Honeywell 7600 recorder, and played back onto a Honeywell 1108 Visi- 
corder plot. See Fig. 10 for a typical Visicorder plot. 

The last piece of apparatus is the photographic equipment. The 
shadowgraph system is simply a point light source from a one microsecond 
duration spark and a camera back (no lens or shutter) placed immediately 
behind the cavity section (see Fig. 9 ). Two types of film were employed, 
Polaroid (3000 ASA) and 70mm (1600 ASA Kodak Recording 2475, on a motor- 
driven film transport). The spark is triggered by an accurate (+ .01 millisec) , 
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adjustable delay circuit. The delay circuit is activated by the initial 
pressure pulse from the chamber transducer. Because the flow under study 
was shown to be reproducible, it was reconstructed by pictures of differ- 
ent test runs taken at progressively longer delays. This method was 
judged superior to a high-speed camera (6000 frames per second) because 
more pictures per cycle and better clarity could be obtained. 

With this apparatus, it is possible to vary many parameters. 

The versatile design of the chamber and cavity make it possible to vary 
apparent shock frequency (defined as the frequency based on the observed 
wave passage time), cavity frequency, orifice diameter, and apparent shock 
strength (pk-pk amplitude at cavity section which can be varied by placing 
the cavity at a pressure node or ant i-node) . Actual shock strength may 
be varied continuously by adjusting pulse gun or chamber pressure. Also, 
more than one gas may be used in the system. 

In this preliminary study, only a certain few parameters were 
varied. The cavity located was in three positions: 

1. at middle of foot-long chamber (i.e. , pressure node) 
with ~ 1100 Hz observed 

2. at end of foot -long chamber (pressure anti-node) 
with ~550 Hz observed 

3. at end of six- inch-long chamber (pressure anti-node) 
with ^1100 Hz observed. 

Shock strengths investigated covered a range of £ values (ratio of pk-pk 
first wave amplitude to steady-state chamber pressure) from .25-1.5, where 
average chamber pressure was held at ~ 68 psia. Initial chamber pressure was 
65 psia for all test runs. The average apparent shock frequency (based 
on the first eight waves) was kept within 10% of the observed natural 
frequency of the cavity. The orifice diameter was not varied. Only the 
nitrogen tests are reported here. 

As has been stated, the data were taken in two forms, pressure 
plots and shadowgraph pictures. The data will be discussed in that order. 

The damping of high amplitude waves by a resonator was studied through 
pressure plots. The effects of two parameters, shock strength and apparent 
shock frequency, were studied. Various shock strengths ( £ = .25, .44, 

.70, 1.0, 1.5) were observed. The results of this study with the cavity 



at the end of the foot-long tube is presented in Fig. 11 . Note that 
there is always a sizable difference between the damped and the 
undamped cases. There is not, however, a clear indication that higher 
amplitude waves are damped more effectively than lower amplitude waves 
or vice versa. Thus, from this series of tests the effect of shock 
strength on resonator damping efficiency does not seem to be dominating. 

To isolate the effect of apparent shock frequency on damping, 
two series of tests were run. In the first series, the cavity (tuned 
for 550 Hz) was at the end of a foot-long chamber with the orifice 
blocked and open. This configuration provided an apparent shock fre- 
quency of 550 Hz. In the second series, the cavity (tuned for 1100 Hz) 
was at the end of a six- inch-long chamber, with the orifice blocked and 
open. In this case the apparent shock frequency was ~-1100 Hz. Note 
that in both series the cavity was at a pressure anti-node and velocity 
node. Equal strength shocks were observed in both series ( £ = .70). 

The higher frequency series (second series) produced very little differ- 
ence between the damped and the undamped cases. Thus, there would appear 
to be a frequency effect influencing the damping efficiency - a point that 
must be checked with a refined shock pulse generator. 

Fig. 12 illustrates a result found early in the study, namely 
that the resonator is more effective at the end of the foot-long chamber 
than at the middle. The "end" represents a pressure anti-node, velocity 
node, and an apparent shock frequency of ~'550 Hz. At the middle, is 
the opposite situation, pressure node, velocity anti-node, and an apparent 
shock frequency of ~1100 Hz. The two preceding paragraphs show that 
a resonator cavity at the end position is more efficient because it oper- 
ates at a lower frequency. Further tests (with the refined shock pulse 
generator) are needed to determine the effect of the velocity in the 
chamber on damping efficiency. 

Aside from pressure data, many pictures were taken of the shock 
wave's interaction with the cavity. To date the most important pictorial 
data provide a detailed survey of a single cycle. The cycle under scrutiny 
is described in Figs, 10 and 13 . Fig. 14 (note orifice outline) shows 
the shock wave entering the field of view. As it proceeds down the chamber 
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to the left, passing over the orifice it develops a "leg" (Fig. 15), 

When the shock is just out of the field of view, a vortex 

ring is visible just beneath the orifice in the cavity. These vortex 

16 

rings are evidence that a jet is streaming into the cavity . Also seen 
in the cavity at this time is what looks like a weak, curved shock wave 
propagating upward from the bottom of the cavity. Perhaps the weak 
cavity shock resulted from the formation of the main (chamber) shock* s 
"leg". Figs. 16 and 17 show vortex rings and cavity shock wave. As 
the shock continues down the chamber, the cavity jet increases in inten- 
sity and the vortex rings become very pronounced (Fig. 18). 

Just after the shock reflects off the wall near the pulse-gun, 
i.e., a little more than half-cycle, the cavity jet dies down. Turbu- 
lent activity is noted at the orifice in the chamber. Figures 19-21 
shows this transition region of the cycle. 

As the shock continues toward the middle, the chamber jet 
intensifies (Fig. 22). Precisely when the chamber jet disappears is not 
clear, but it is not visible when the shock wave is again within view. 

It is noteworthy that the cavity jet seems to last longer than 
the chamber jet. 

Fig. 23 shows the shock (traveling to the right) entering field 
of view. There is what may be the last vestiges of chamber jet activity 
around the orifice. Figure 24 shows the shock a bit later, over the 
orifice. Compare the picture with Fig. 25 which was taken at the same 
corresponding time, with the orifice closed. This comparison leaves 
little doubt that the orifice is responsible for the formation of the 
"leg" on the shock wave. 

This has been a report on the preliminary study of the inter- 
action of high amplitude (shock) waves with a properly tuned Helmholtz 
resonator. Further study will utilize a more sophisticated shock wave 
generator to eliminate this double-wave complexities inherent in the 
present system. This will make possible a quantitative study of the 
interaction. Future work will study the effects of all variable parameters, 
including the use of two gases to more closely duplicate the conditions 
in a hot- flow combustion chamber. 
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IV. ENERGY SOURCES 

A. USE OF SHOCK WAVES TO STUDY COMBUSTION PARAMETERS 

The purpose of this research, which was described in NASA 
CR 72270 1 and in NASA CR 72371, is to investigate the functional forms 
of the mass, momentum and energy sources in a combustion chamber. Thus 
the components in the liquid phase are visualized as sources (or sinks) 
for the components in the gaseous phase. 

The usual way of performing such an investigation would be to 
solve a set of equations and then to compare the solution with experi- 
mental results. To follow this approach it is necessary to have as many 
equations as unknowns. Since the number of basic equations are always 
less than the number of unknowns, extra equations would be needed. These 
added equations would reflect assumed models for the sources to be in- 
vestigated. Furthermore, a large number of equations would finally 
result and additional simplifying assumptions would generally be needed 
if they are to be solved. The final comparison between analytical and 
experimental results would not be too conclusive, particularly so if 
discrepancies between analytical and experimental results are found. 

Another way of performing this investigation is to use directly 
the experimental results to determine some of the unknowns instead of 
adding extra equations. Thus one ends up with a smaller number of equa- 
tions which satisfy "exactly" the experimental results. Moreover, having 
a smaller number of equations one can usually solve them more accurately. 
Solutions for the functional forms of the sources are found in this 
approach, instead of assumed and then verified. Thus in the direct use 
of measured parameters, one exchanges known information about these param- 
eters for unknown information about the sources. This direct approach 
therefore requires a smaller number of assumptions , leads to fewer and 
simpler equations, and is more accurate. 

We chose to follow this second, direct approach but unfortunately 
there are three sources (of mass, momentum and energy) and it is difficult 
to measure directly and completely three parameters. Thus it was necessary 
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to settle for a compromised approach. It may be noted that, in the 
direct approach, the accuracy and reliability of the measurements are 
very important, but, in principle, they are no more important than in 
the indirect approach. 

Actually instead of directly measuring some combustion param- 
eters, a shock can be induced in the combustion chamber and some of 
the shock parameters can be measured. Information on the combustion 
parameters can be gained by studying the propagation into the combustion 
chamber of this externally generated shock. Consider a long combustion 
chamber where one -dimensional steady-state, distributed combustion is 
assumed to occur. The parameters of a shock, started, say, at the 
nozzle end, will depend on the steady-state condition of the undisturbed 
gas ahead of the shock as well as on the changes of the combustion 
process which occur after the shock front. These changes are primarily 
due to the particle velocity induced by the shock if the shock is of 
moderate strength and droplet shattering is excluded. 

More specifically, through the shock front the Rankine-Hugoniot 
equations hold regardless of whatever mass, momentum and energy additions 
(or subtractions) occur after the shock front. Thus the use of the 
Rankine-Hugoniot equations through the shock front together with the 
steady-state conservation equations for the undisturbed gas between the 
injector end and the shock front, can be applied to yield information 
on the steady-state combustion parameters. This part of the problem is 
here referred to, as the study of the "steady-state" combustion. On the 
other hand, suppose now that the steady-state parameters of the gas 
ahead of the shock front are known, together with the parameters of the 
shock front, and of the shocked gas at some initial time, then the values 
of these parameters could be calculated at any later time, if the 
functional forms of the mass, momentum and energy sources were known. 
Conversely, if three of these parameters were measured in their time and 
space domains, we could deduce the functional forms of the mass, momentum 
and energy sources. This part of the problem is here referred to, as 
the study of the "unsteady-state" combustion. We will measure only pres- 
sure and shock velocities, thus some extra assumptions will have to be 



61 


made. However, the goal is still that od deducing the functional forms 
of the sources based upon some directly measured parameters. Notice 
that before attempting solutions to the "unsteady" part of the problem 
it is necessary to solve for the "steady" part first. 

The results of the steady-state calculations will presently 
be discussed together with some considerations on the equations for the 
unsteady part. But first a brief description of the experimental set up. 

The Experimental Apparatus 

A liquid propellant rocket motor with square cross section 
(area equivalent to a circular section of 3-inch diameter) is the prin- 
cipal piece of test apparatus in this research. The motor is made up of 
sections which can be added so as to form a chamber up to more than 20 
diameters in length. Transducers are then inserted at appropriate inter- 
vals with as many as six Kistler piezoelectric crystal transducers used 
simultaneously to record local pressures. Rocket tests are expected to 
be reproducible and they are checked for this property via repeated 
testing. 

Further details of the design are as follows: the injector 

utilizes a 16-doublet ,like-on-like design to produce axially-distributed 
combustion which is an important aspect of the investigation. The shock 
itself is produced by a shock tube mounted at the nozzle end of the cham- 
ber, facing the injector, with the nozzle consisting of 12 individual 
radially oriented nozzles (see Fig. 26). The diameter of the shock 
tube is equal to the diameter of the motor so as to facilitate the gen- 
eration of a simple planer shock front. Nitrogen is normally used as 
the pressurant in the shock tube. If the combustion chamber were filled 
with inert gases the shock generated by the shock tube would be an almost 
perfect step function moving at constant speed, however, with the combus- 
tion present the shape and velocity of the shock is altered. Through 
the study of these alterations we are attempting to infer the functional 
form of the sources due to the presence of the liquid propellants and of 
the combustion. Undoubtedly, the production and recording of the data 
which is required, with the accuracy which is necessary, is a difficult 

* Pressures to 1500 psig have been used and together with thin stainless 
steel burst disks (.008 to .015 inches) the desired shocks are generated. 
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Figure 26 


Shock tube combustion chamber for studies of combustion response 
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experimental task. Up to now the success has been limited, but 
encouraging. 

The measurement and recording techniques, which have been 
necessary to obtain suitable data, are representative of the current 
state of the art. Pressure measurements require flush-mounted pickups 
with high frequency response since the waves to be measured are 
steep- fronted. The most successful approach has been the use of 
ablative coatings to protect the diaphragms of Kistler 603A transducers 
which possess a natural frequency in excess of 400 kHz. 

Protection of Kistler-type piezoelectric transducers was first 
demonstrated by Rogero at JPL^’ The tests described here use 
similar arrangements with the addition of a water cooled-housing 
to prevent heat soak and to limit the thermal drift during short dura- 
tion rocket firings and during any prerun temperature variations. The 
water passage simply surrounds the transducer to eliminate heat transfer 
paths other than through the ablative itself. 

Side-mounting of the transducer in the square-motor allowed 
typical shock transient times greater than ten microseconds, which was 
a less stringent condition than the transient times of five microseconds 
or less on the wall of an actual rocket chamber. On -the -other -hand, 
the shock wave presented a waveform that was steeper than many of the 
wave shapes encountered in an unstable rocket. Although not a necessary 
condition for these tests, end-mounting presented the limiting situation 
of zero transient time (all portions of transducer face received the 
pressure signal simultaneously) and the wave was a steep-fronted shock 
of high strength. Based on these differences it was found in the tests 
that side-mounting of the transducers was far less critical than end- 
mounting from the standpoint of initial overshoot and cavity resonance. 

By cavity resonance, reference is to the cavity separating the transducer 
diaphragm from the chamber wall which has an associated frequency spectrum. 
The cavity resonance was evident in only the tests without ablative and 
using the .125 inch recess (this ablative thickness is routinely used 
on the transverse mode tests where the wave shape is often less critical 
from a measurement standpoint and the heat transfer environment is more 
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severe. The side-mounted, thick ablative coating still exhibited 
an initial overshoot of 20%, however, while the thin ablative and 
flush-mounted configurations experienced no overshoot when similarly 
mounted . 

The ringing frequency of the transducer following the shock 
arrival was used as a measure of the natural frequency of each unit. 

The flush -mounted transducers provided the 400 kHz basis for compari- 
son. The thin ablative coating indicated frequencies of 300 kHz while 
for the case of the thick ablative the transducer natural frequency 
was reduced to 90 kHz. These data are summarized in Table I. 

The point to be made here is that even the thick ablative 
approach will be found to often possess better frequency character- 
istics than many typical transducer types and mountings. A 90 kHz 
natural frequency is still quite high when compared to many strain 
gage and capacitance- type transducers which are flush-mounted but 
require water-cooled diaphragms. Therefore in limited duration firings 
(less than five seconds) the use of suitable ablative coatings may prove 
to be highly beneficial - especially when analyzing higher modes, and 
nonlinear or shock-type waves. 

One other transducer approach that has proven successful to 

18 

researchers at the University of Michigan has involved the use of a 
ferroelectric element. Specifically, the use of lead metaniobate (with 
a tin backing rod designed for removal of the reflected wave within the 
sensor) is the basis for the transducer design. This approach is being 
carried along as an alternate. Testing has revealed the design to be 
less rugged than the Kistler transducers although the initial wave shape 
recording is quite good. Later testing may still require this approach. 

Returning to the other principal factor in successful data 
collection, the recording system, the approach utilizes a Honeywell 7600 
tape recorder. This equipment limits dynamic time displacement error 
between channels to less than jh 1 microsecond. Through the use of the 
optimum transient response filters in the 1 7/8 inch/sec playback, over- 
shoot is reduced to a minimum level (order of 5%). The function of the 
double extended response mode allows an 80 kHz band width at a tape speed 
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of 120 inches/sec. and a small sacrifice in the signal- to-noise ratio 
(order of 3 db) , Signal conditioning equipment of the charge amplifier 
type was specifically designed for these tests. The equipment is 
capable of conditioning the signals from a wider range of transducers 
than standard charge amplifiers (e. g. , both the piezoelectric and 
ferroelectric transducer types can be accommodated). The history of 
the overall recording arrangement has been one of small improvements 
at a number of points in the system resulting in the final desired data. 

Fig. 27 shows the pressure traces associated with the helium- 
filled chamber tests previously mentioned. Both time of arrival data 
and step amplitude are recorded to- the desired accuracy. For comparison 
the same chamber with combustion taking place presents a somewhat differ- 
ent situation as shown in Fig. 28 . The downstream pressure traces 
are very similar to the inert gas tests. However, as the shock wave 
moves into the region of greatest interest, that zone in which the major 
portion of the combustion is taking place, the step-like shape is con- 
siderably altered. Pressure is seen to rise following the passage of 
the wave and there are considerable pressure variations present. These 
variations are not because of transducer ringing or other instrumentation 
effects but are a function of location, and hence a function of the com- 
bustion conditions present. 

It is from records such as these, together with the theory 
presented, that the alterations from the ideal case will be analyzed and 
interpreted to explain the steady-state and finally the unsteady combus- 
tion phenomena. 

Calculation of the Steady-State Combustion Parameters 

In last year's progress report"^ (NASA CR 72270), equations were 
given which allow the calculation of the steady-state combustion parameters 
after having measured the static pressure in the combustion chamber and 
the velocity and overpressure of an externally induced shock wave. For 
the derivation of those equations it was assumed: 

1) steady state between the injector end and the station of 
interest 


2 ) 


one-dimensional flow 






- 300 psig 
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3) instantaneous vaporization of the oxidizer and distrib- 
uted vaporization of the fuel 

4) no change in the mass, momentum and energy of the drops 
through the shock front itself 

5) plane shock front 

6) boiling temperature of the fuel drops 

7) instantaneous reaction of gaseous fuel and oxidizer to 
give equilibrium reaction products. 


Notice that no droplet drag and evaporation model was assumed. The measured 

/ 

shock velocity and shock overpressure complete the system of equations, 
otherwise, two more relationships would have been needed. 

Those equations are for a fuel-oxidizer system made up of 
C, H, 0 and the products are allowed to contain the following eight 
compounds: OH , CO , C0£ , 0 Oj , H ^ > Hg-O , H . Since the eight 
combustion products are assumed to be in equilibrium before and after 
the shock front we have 16 chemical equilibrium equations, which added 
to the three conservation equations from the injector to the shock front, 
and to the three Rankine-Hugoniot relations through the shock front, 
give a total of 22 nonlinear algebraic equations to be solved simulta^ 
neously. The solution of the system presented some difficulty but 
finally it was achieved, however, not by the method indicated in last 
year’s progress report. The actual solution was accomplished by using, 
as a parameter, the factor ^ A which gives the mass flow per unit 
time of the liquid drops at any given station ( f = mass of liquid per 
unit volume of the combustion chamber, u^ = average liquid velocity, 

A = cross section area of the combustion chamber). In this approach the 
shock velocity (U) is among the calculated parameters. We can then plot 
f u^ A versus U for any given pressure ratio through the shock front . 

The measured U will then supply information as to the actual value of 
p u, A thus determining the entire solution. Actually the solution of 

h 1 

the system is very sensitive to the value of ^ u^ A . Only within a narrow 
range of j^u^ A is the velocity of liquid less than the velocity of 
the gas and greater than zero. Hence the solution of the system is 
determined within an acceptable range even without the measurement of U . 

The measured U should then fall within the calculated range. The 
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measured static pressure along the combustion chamber was also used 
directly to obtain the above solution. The reason for the leading role 
played by u^ A is its controlling effect in the energy equation. 

The amount of fuel in the liquid phase does not greatly influence the 
momentum equation. It influences the mass equation. But, most of all, 
it influences the energy equation since it control Is the amount of 
chemical energy released and this, in turn, is what makes the entire 
engine work. 

The results may be discussed with the help of Figs. 29 and 30. 


These data were calculated for ethyl alcohol and liquid oxygen; fuel 
mass flow = 428. g/sec • oxidizer mass flow = 986 g/sec , and static 
pressure in the engine at the injector = 296. psi. In Fig. 29 the 



distance from the injector. In the table the values of the gas velocity 
(u^) and liquid velocity (u^ ) for each value of the above fraction are 
noted. For instance, at 30" from the injector the acceptable value of 
A j^u ? / Wp ranges roughly between .29 and .23 . For .29, u = 10,090 cm/ sec 
and u^ = 14,030 cm/ sec ; for .23, u =45 cm/ sec and u^ = 14,730 cm/ sec . 
It can be concluded, then, that A q u t / cannot be much less than .23, 
otherwise the liquid would be moving toward the injector. This term 
cannot be much more than .29, otherwise the liquid would be moving faster 
than the gas at a distance of 10 chamber diameters. Similarly, the 
other acceptable values of the combustion parameters vary in narrow ranges. 

We have already seen that u^ must be roughly between 14,030 and 14,730 
cm/sec . From Fig. 29 we see that the temperature (T) must be between 
3199°K and 3236°K . The selection of the exact value of A ^ u fc / ui p is 
to be made by the measured shock velocity for a given shock front over- 
pressure. From a practical viewpoint, the values of the acceptable shock 
velocity also varies in a narrow range (within the limit of the experi- 
mental error) and therefore it is then only necessary to verify that the 
measured velocity is compatible with the calculated one. In Fig. 30 we 
have the measured static pressure loss which was used directly in the 
calculations, and the approximate acceptable ranges for the values of 
the gas velocity (u^) , the gas density (j^) and the chemical energy 
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released per unit weight of the combustion products (H^) , The dashed 
lines within the ranges are the likely values under the assumption 
that the velocity of the liquid drops approach the gas velocity as 
the distance from the injector increases. 

The most important of the findings is that one (or more) of 
the stated assumptions is not valid for the region between the injector 
and about 3 or 4 diameters from it. It was assumed that the oxidizer 
vaporizes immediately and completely (Assumption 3) and that instantaneous 
reaction of the gaseous fuel to equilibrium combustion products occurs 
(Assumption 7). Thus, very close to the injector, a very low calculated 
temperature was expected, since there the fuel to oxidizer ratio is 
very unfavorable. We expected this region to extend out to, perhaps, 
one diameter from the injector and we were ready to accept the failure 
of the model in this region. In other words, we expected the temperature 
curve of Fig. to be shifted to the left by approximately two diameters. 
The actual findings are in contradiction with expected ones for a dis- 
tance from the injector too long to be acceptable. A critical review 
of the 7 assumptions made, in the light of the theoretical and experi- 
mental results indicates that Assumption 1 is verified, Assumptions 4, 

5 and 6 are not very influential so that one or more of Assumptions 2, 

3 and 7 should be responsible for the calculated results. It is felt 
that the assumption of immediate vaporization of the oxidizer (Assumption 
3) is the weaker one. Since some fuel is calculated to exist in liquid 
form even at 10 diameters from the injector, it is not unreasonable to 
expect some oxidizer to exist in liquid form out to 3 or 4 diameters 
from the injector. 

Accordingly, we have modified our equations to account for the 
presence of liquid oxidizer in the vicinity of the injector. We have 
assumed that the average velocity of the liquid oxidizer is equal to that 
of the liquid fuel and we have added a new variable, expressing the mass 
of liquid oxidizer per unit volume of the combustion chamber. Then one 
more quantity needs to be measured. We are now using two shocks, one 
from the nozzle and one from the injector. Measurements of these shock 
velocities and overpressures are sufficient to obtain solutions to this 
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case involving liquid oxygen near the injector. However, it turns out 
that it is not convenient to solve the new set of equations for the 
expected values of the experimental results and then to select the actual 
solution later (as it was done for the no-liquid-oxidizer case) . We 
must wait for sufficient experimental results before performing this 
second calculation. The solutions given in Figs. 29 and 30 are still 
expected to hold for distances from the injector greater than 3 or 4 
diameters. 

Other interesting steady-state results are: 


1) If the static pressure drop from the injector to any station 
is taken to be of order 1, then the change of momentum of 
the gas between the injector and that station is of order 1 
and the change of momentum of the liquid is of order 10 ^ 
(combined effect of drag and vaporization) . In other 
words, for the calculation of the steady state the influ- 
ence of the liquid in the momentum equation of the gas 

is negligible: changes of momentum of the gas are balanced 

almost exclusively by changes in static pressure. 

2) The factor A ^ u^ is determinable due to its strong 
influence on the energy equation while f and u^ cannot 
be determined due to their weak influence on the momentum 


3) 

4) 

5) 


equation. 

The energy source is a function of the mass source of the 
local fuel to oxidizer ratio and of the local temperature. 
Local fuel to oxidizer ratio and local temperature are then 
related to each other through the hypothesis of chemical 
equilibrium. 

Y (ratio of specific heats) varies along the chamber 
approximately from 1.27 (near injector) to 1.22 (far from 


injector) . 

If the local speed of sound is calculated with 



even using the local y , the error might be as high as 


10% due to having neglected the chemical processes. 
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6) The chemical energy released (H^) is approximately 
redistributed as follows: 

65% internal energy (warming up of the combustion 
products from the liquid temperature to the 
combustion chamber temperature) 

20% mechanical work (P/f) 

7% vaporization of fuel 

7% vaporization of oxidizer 

kinetic energies of the gas and of the liquid 
are negligible 

Summary 

In conclusion, two sets of equations were developed for the 
computation of the steady-state combustion parameters using directly 
measured values of overpressure and velocity of an externally induced 
shock wave. One set of equations is valid for regions where only 
liquid fuel is expected to be present. This set of equations can 

actually be solved without experimental data and such data can then be 

used a posteriori. A specific set of solutions are presented in Figs. 

29 and 30. Although complete experimental data have not yet been 
collected, preliminary experimental data fall in the expected calculated 
range. The second set of equations is valid for regions where both 

liquid fuel and liquid oxidizer are expected to be present. For the 

use of this second set, experimental data must be available a priori. 

We are presently gathering these experimental data. However, this second 
set of equations has already been solved using assumed values for the 
variable to be measured. To complete the steady-state solution it will 
be necessary to gather sufficient experimental results, and these tests 
are presently in progress. 

The study of the unsteady-state has also been initiated. 

Great effort has been concentrated in the selection of the proper set 
of assumptions and of equations in the light of the preliminary results 
of the steady-state portion and of experimental data. It should be evi- 
dent that final solutions are meaningful only if consistent assumptions 
and proper equations are used. We are presently nearing the end of this 
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selection process. The next step will be the solution of the selected 
set (or sets) of equations. Simultaneously, experimental data for the 
unsteady part of the problem will continue to be gathered. Perhaps the 
most severe experimental problem encountered in recent weeks has been 
the occurrence of spontaneous instability becaused of the reduced 
damping associated with the radial nozzle, shock tube design. The 
limited use of acoustic dampers at the nozzle section has already been 
shown to be highly successful in eliminating this problem. 


B. CHEMICAL KINETICS AS A DRIVING MECHANISM FOR LIQUID PROPELLANT 
ROCKET COMBUSTION INSTABILITY 

The Combustion Process and Combustion Instability 

It is widely accepted that high-frequency combustion instability 
consists of organized oscillations of the properties within the chamber, 
which are maintained and amplified by the combustion process itself. 

The combustion chamber gases are put into a forced oscillation driven by 
the combustion process interacting with the resonance effects of the cham- 
ber geometry. The resulting oscillation patterns are very similar to 
those of the acoustic modes of the chamber. Thus, there seem to be 
coupling mechanisms between the combustion process and the acoustic prop- 
erties of the chamber which can organize and amplify random perturbations 
into periodic oscillations and then maintain these organized oscillations 
at a significant amplitude. Questions arise as to the nature of these 
coupling mechanisms involving the combustion process: how they organize 

and drive the oscillations, and which phase of the complex combustion 
process they are associated with. 

The "combustion process" in a liquid propellant rocket engine 
involves everything that happens to the propellants between the instant 

they are injected into the chamber and the moment conversion to gaseous 

/ 

products is completed. Involved are atomization, vaporization, mixing, 
burning, recirculation, and many other steps. Since it is the chemical 
reaction of burning (fast oxidation) that supplies all of the energy 
generated within the chamber, it seems logical that this process also 
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supplies the energy for the organized oscillations of instability. 

The problem is how this energy release is related to the organized 
acoustic-mode oscillations. 

The rate of energy release depends on the rate of chemical 
reaction. For example, if the rate of the reaction - the kinetics - 
varied periodically, so would the rate of energy release. Such effects, 

involving when and where energy is released in the chamber, could have direct 
bearing on the nature of the oscillations produced. It was felt that 
further study of the nature of a kinetic coupling mechanism should be 
conducted to provide a fuller understanding of combustion instability 
and of the combustion process itself. 

Joint Study 

During the past year, a program has been implemented at 
Princeton to investigate such an assumed coupling. This program attempts 
to determine what relationship, if any, chemical kinetics has with com- 
bustion instability associated with liquid propellants, and the 
nature of that relationship with regard to factors such as phasing and 
amplitude. 

Because this investigation would study both the combustion 
process and combustion instability, it was implemented by both the 
Combustion Group, under Professor Glassman, and the Liquid Rocket Group, 
under Professor Crocco and Mr. Harr j e , as a joint project. The essence 
of the investigation, then, is to combine both avenues of research, com- 
bustion processes and liquid rocket instability, to glean a broader know- 
ledge in areas of common interest. 

The Combustion Group has had much experience in pure kinetic 
studies. One important piece of equipment developed by that group is 
the flow reactor which is being used in fundamental studies of oxidation 
(combustion) reactions. The reactor spreads out the oxidation reaction 
in both time and space, thus allowing precise measurements of kinetic 
data. Hydrazine and nitrogen tetr oxide oxidations have been studied 
extensively, and hydrocarbon studies are just beginning. Also the 
Group is studying combustion instability in a gas rocket, where all liquid 



effects and changes of state are nonexistent. Thus the combustion 
process is much simpler, and the kinetics are exposed for study without 
the added complexities of atomization, vaporization, etc. The Liquid 
Propellant Rocket Group provides rocket test facilities and long expe- 
rience in the study of combustion instability in liquid propellant 
engines. It was felt that a combination of these three methods of 
investigation in order of increasing complexity - reaction kinetics, 
gas rocket instability, and liquid rocket instability - would carry the 
study step-by-step to a useful conclusion. What would be learned from 
each step could be applied to the next one to advance knowledge in a 
systematic manner. A unified treatment could be achieved because the 
same chemical systems could be studied from each point of view, and 
data from one approach would be directly applicable to the others. 

Methane-Oxygen Propellant System 

The investigation has begun with the simplest hydrocarbon 
propellant system: methane (CH^) and oxygen (Op . In addition to the 

propellant simplicity, this propellant combination is related to the 
hydrocarbon fuels previously tested in the instability research, and 
possesses performance characteristics that make it an attractive choice 
for future propulsion applications. 

Concerning the interrelated studies, the flow reactor has just 

undergone extensive modifications; its first investigation will be the 

kinetics of methane oxidation. (The modifications were necessary because 

of the unusually long induction time - of methane compared with other 

hydrocarbons and other compounds, such as alcohols.) Some experimental 

work has been done with gas rocket instability of methane-oxygen propel- 
19 

lants and more is planned. Thus the kinetics and gaseous combustion 
of methane should be further clarified within the near future. The 
liquid propellant rocket facilities have been modified so that methane 
may be used as a propellant. The me thane -oxygen system in the liquid 
program becomes, of course, liquid methane - liquid oxygen. We shall 
abbreviate this propellant combination as LOX/LCH^. Both substances in 
their liquid forms are cryogenics. The normal boiling point of LCH^ is 


* the delay between injection and reaction initiation 
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- 259°F (201°R) . Its freezing point is - 298°F (162°R) , so its liquid 
range is small. This causes a number of difficulties in its use. 

From the particular point of view of the liquid rocket test 
program, the LOX/LCH^ propellant system offers an additional check of 
present theoretical models, and a further tool for unraveling experi- 
mental mysteries. The double-cryogenic propellant combination offers 

a new system to study from an n - H standpoint, and further investi- 

11,20 

gation of the velocity effects model of transverse instability. 

Also, this system provides a firm base on which to begin kinetics studies 
in the rocket motor and to investigate the puzzle of the direction of 
spin of transverse modes of instability. 

Spin Direction Dependency Upon Fuel 

It was decided that the point of departure of the experimental 

program would be an interesting series of previous tests concerned with 

the spinning first tangential (IT) mode of instability: it had been 

found that for two different propellant combinations, liquid oxygen and 

RP-1 , and liquid oxygen and ethyl alcohol, the directions of spin of the 

IT mode were different. The spin directions were in opposite directions 

for the two propellant systems; that is, with the same injector -chamber 

configuration, one system had instability spinning clockwise, while the 

3 

other system had instability spinning counterclockwise . These results 
were repeatable, and evidently they represented characteristics of the 
propellant systems. Here, then, was an effect that clearly might be 
related to the combustion kinetics. Therefore, the investigation was 
begun at this point in an effort to discover the reason for such an effect 
due to propellants . 

Nomenclature 

It will be helpful to introduce some nomenclature at this point. 
Liquid oxygen is abbreviated LOX, and ethyl alcohol is abbreviated ALC. 

F stands for fuel; 0 stands for oxidizer. The injector used the "spud" 
design with each spud containing injection orifices. The spuds could be 
like-on-like impinging (L-L) or unlike impinging (U) . An L-L spud has 
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two pairs of holes, one pair for fuel, one pair for oxidizer. The 
spray fans formed by the impinging liquid streams overlap. The line 
of centers of the orifices may be oriented tangentially (T) , or 
radially (R) , with respect to the chamber geometry. The fuel orifices 
are placed nearer the wall because of heat transfer considerations. 

The pairs of orifices must be displaced because of manifolding require- 
ments. A typical L-L, T spud is shown below 


^Wall 

JJJJ-Lkdjji 



A bar is placed over the F to indicate that the fuel orifices were 
nearer the wall: F . The direction of spin of the instability is noted 
by the order in which the pressure wave crosses the respective orifices: 

0 — > F (clockwise) or F 0 (counterclockwise) . 

Previous Testing and Planned Program 

The results of the previous tests on spin direction had been 

3 

the following: LOX/ALC spun 0 — ¥ F; LOX/RP-1 spun F — * 0 . It was 
reasonable to presume that kinetics could have something to do with the 
spin direction because ethanol and RP-1 are fundamentally different types 
of molecule. Ethanol is an alcohol, with the specific functional group 
hydroxyl; RP-1 is a mixture of hydrocarbons with no functional groups. 

But not only simple chemical differences exist between these two sub- 
stances. There is a considerable difference in molecular structural com- 
plexity, and therefore presumably, other than chemical differences exist 
between them. These differences could be minimized by using a hydrocarbon 
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and its corresponding alcohol, methane and methanol, ethane and ethanol, 
for example, as fuels. Thus, a better assessment of the importance 
of the functional group could be made. Additional tests could then be 
made with other fuels, with chemical properties changed step-by-step 
by chain lengthening and side chain addition. Hopefully, this series 
of tests would isolate the chemical effects and their influence on 
instability by correlating changes in instability pattern with changes 
in fuel system. Clarification of the coupling mechanism would be further 
enhanced if variations in physical properties associated with the combus- 
tion process; e.g., surface tension, heat capacity, boiling point, and 
heat of vaporization, were minimized. 

Liquid Methane Facility Modification 

Since the first series of liquid rocket tests were made using 
the LOX/LCH^ propellant system, the first requirement was the modification 
of the liquid propellant rocket test facilities for use of liquid methane 
fuel. Liquid methane is cryogenic, and while the facility had for many 
years used cryogenic liquid oxygen, a cryogenic fuel presented additional 
problems. The boiloff gases are combustible with the ambient air, and 
their low density allows them to collect in pockets near ceilings if care 
is not taken in the same manner as with liquid hydrogen. 

The liquefaction of the methane and its storage in liquid form 
presented problems also. To buy liquid methane would not only be expensive 
but also inconvenient from a delivery and storage standpoint. It was 
decided to buy gaseous methane and then liquefy it at the facility. 

This decision required the design and construction of a methane condenser 
with adequate capacity for the projected use rate of methane. A simple 
coiled tube, through which gaseous methane is flowed, is immersed in a 
bath of liquid nitrogen. This arrangement provides ample liquid methane, 
and at a factor of four savings over commercial purchase. Even so, the 
liquid methane costs about $2.50 per pound, so it cannot be wasted. A 
vacuum- jacketed pressure tank is used for runs, and it is able to store 
liquid methane for sufficient periods of time. The tank temperature is 
controlled by flowing liquid nitrogen into the vacuum jacket and then 
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draining when the tank is cold enough. Because the liquid methane is 
so expensive, liquid nitrogen is used for all other cooldown opera- 
tions as well. Since liquid nitrogen is cold enough to have the capa- 
bility of freezing liquid methane, care must be taken during the con- 
densation process that the methane flow rate does not become too low. 

It appears that all problems with the liquefaction process, cooldown, 
and storage have now been solved. Routine procedures have been devel- 
oped for every phase of the operation. 

Results of First LOX/LCH^ Tests 

So far, thirteen runs have been made with the LOX/LCH^ pro- 
pellant system. Four injector configurations have been used, and a 
three-blade baffle was added to three of them. Mixture ratios varied 
between about 1.9 and 4.0; this range corresponds to a range of equiv- 
alence ratios (r/r ..... . . ) between 0.5 and 1.0. Chamber pressure 

stoichiometric 

was nominally 150 psia. 

All tests showed spontaneous instability growing out of the 
combustion noise. The oscillations began randomly, with essentially 
zero magnitude and organized into periodic oscillations while the 
amplitude increased to significant values. The growth of a typical 
oscillation can be schematically pictured as 



The oscillations always began smoothly; there never appeared any pulse 
or spike which might have triggered them. 

The mode was mostly the first tangential (IT) spinning mode, • 
although four runs showed basically 2T modes. These four runs employed 
a three-bladed baffle. The spin direction was always F —►0 . This spin 
direction characteristic of LOX/LCH^ is the same as the direction pre- 
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ferred by LOX/RP-1, and opposite to that of LOX/ALC. Note that both 
RP-1 and methane are hydrocarbons. Thus their similar behavior, as 
contrasted with that of alcohol, is not surprising. 

Four runs showed mixed modes, with basically 2T mode, and 
frequencies higher than the 2T superimposed over it. The instability 
damped out of its own accord in four other runs in which the three- 
blade baffle was used. In one run, the motor was pulsed in the 0— *F 
direction; however, after the pulse decayed, the instability returned 
to its original direction, F-*0, with about the same amplitude as 
before. 

Peak-to-peak amplitudes ranged between 20 psi and 420 psi, 
the magnitude depending on the injector system and the presence of the 
baffle. There appears to be no correlation between either amplitude 
or frequency and mixture ratio; rather, both appear to correlate with the 
injector-baffle system. 

The first four runs, 2227-2230 were made with a 1 x 16 like- 
on-like impinging doublet injector in the tangential orientation. In 
the past, this injector has exhibited stable operation with the ability 
to produce nonlinear instability with other propellant systems. All 
instability in these tests was spontaneous, in the IT spinning mode in 
the F 0 direction. The instability did not damp out. The correspond- 
ing amplitudes and frequencies are detailed in the table to follow. 

These amplitudes and frequencies are given for well -developed instability: 
as the instability grows out of the noise, its amplitude and frequency 
normally increase together until they reach roughly fixed values. 

The next four runs, 2238-2241, were made with a 6 x 2 mixed 
spud injector, alternating groups of two spuds with like-on-like impinging 
doublets with groups of two spuds with unlike impinging doublets. A 
three-bladed baffle was also used. Again, instability was spontaneous, 
and it did not damp out. It was of the spinning type, F-^0, but the 
mode was mixed. The mode was mostly second tangential, but higher fre- 
quency modes were superimposed over that, often obscuring it. A pattern 
of development of these higher frequencies was noticed. When the mode 
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was primarily 2T, at the beginning of the run, it was usually of higher 
frequency, and when the higher overtones became prominent, the 2T 
frequency decreased. Amplitudes varied with the frequency, as shown 
below. 

The following three runs, 2242-2244, were made, again with 
the three-bladed baffle, but with a 1 x 12 like-on like impinging 
doublet injector. This injector has shown more stable operation than 
the 1 x 16 injector in the past. Instability was still spontaneous, but 
in these runs, it damped out by itself before the end of the run, and 
the remainder of the run was stable. The mode was IT, spinning, F 0 . 

The frequencies were all around 2500 cps, and the maximum amplitudes 
were all around 40 psi. Damping took 0.34 to 0.53 second after reaching 
maximum amplitude. The damping times were about ten times as long as the 
times taken for the instability to develop fully. 

The last two runs, 2245 and 2246, were made with the 1 x 12 
like-on-like injector from the annular motor. This injector uses 
recessed doublets with a smaller impingement angle'*' and is extremely 
stable, with associated lower performance. The three-bladed baffle was 
again used. The instability was spontaneous in the IT, spinning F-^0 
mode, but it attained an amplitude of only about 25 psi and a frequency 
of around 2000 cps. Run 2245 damped out very quickly, in about 70 millisec. 
after reaching full development. Run 2246 was pulsed with the gun in the 
0 -*>F direction. Before the pulse, it had been spontaneously unstable 
in the F-^0 direction. After the pulse, it resumed its F 0 direction 
at about the same amplitude. It seemed to begin damping toward the end 
of the run. 

The results of these tests are conveniently summarized in the 
following table. Question marks in the JT column indicate that the mixture 
ratio was unknown due to a malfunction of the metering system, which could 
not be repaired once a series of tests was in progress. The value entered 
is the desired quantity. The true value was probably fairly close. 

Question marks in the Frequency column indicate that analysis was diffi- 
cult because of the mixed modes. 




LOX/LCH PROPELLANT SYSTEM, 9 -inch MOTOR, RUN DURATION 1 second 



























85 - 


Analysis and Conclusions 

The liquid oxygen-liquid methane propellant system obviously 
has a high tendency toward instability. Every test so far has shown 
linear instability spontaneously growing from combustion noise into 
organized oscillations. In all tests, the transverse instability 
appeared in the spinning form, and the direction of spin was always 
F—+ 0 . This spin direction is evidently characteristic of LOX/LCH^, 
and strongly preferred by it, as shown by the inability of a pulse 
from the opposite direction to reverse the spin. The F-*0 tendency 
overcame the applied pulse direction to return to its original direction. 
The mode, frequency, amplitude, and damping observed seemed to be depen- 
dent upon the injector used and the presence of the baffle. Thus the 
characteristics of the propellant system alone are taken to be the 
spontaneity of instability, its appearance in spinning form, and the 
direction of this spin. These are the factors which may be compared 
with those of other propellant systems to determine which part the combus- 
tion process plays in driving instability. 

A comparison may be made between these results and previous 
results obtained in the same motor with two other propellant systems, 
LOX/RP-1, and LOX/ALC. All three series of tests used the tangential 
spud orientation. Hence, all the systems were sensitive to the mixing 
effects of the tangential velocity associated with the spinning wave.'*"'*' 
The tangential velocity effects are always present in spinning modes. 

The velocity effects referred to here involve movement of one propellant 
relative to the other, thus setting up regions of unequal mixture ratio 
which therefore have different energy release rates. That propellant 
which vaporizes more easily is more likely to be moved relative to the 
other, for the gaseous component of the propellants will respond most 
readily to pressure and velocity perturbations. This effect is inten- 
sified if the vaporization rate of the second propellant is so slow that 
its droplets are still very large and thus are not carried along with the 
gases. It is possible that this effect helps to determine the spin 
direction. Accordingly, vaporization data for LOX, RP-1, ethanol, and 
liquid methane are presented on the next page. 
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PROPERTY 

LOX 

ETHANOL 

RP-1 

LCH. 

4 

Normal Boiling 
Point , °K 

90 

351 

490 

111 

Thermal Conduc- 





... ... Btu 

tivity 

hr ft °F 

0.08 

0.094 

0.077 

0.11 

Heat Capacity 

0.41 

0.60 

0.5 

0.85 

Btu/lbm °R 





Heat of 

Btu 

Vaporization, lbm 

92 

362 

125 

219 

Spin Direction 

— 

0 -►F 

F — * 0 

F-*" 0 


Notice the low boiling point of liquid oxygen, and its small 
heat capacity. Therefore, it will reach its boiling point very quickly, 
and its low heat of vaporization assures that it also will vaporize 
quickly. For these reasons, one can confidently predict that in any 
propellant system involving LOX and one of the listed fuels, the LOX will 
vaporize first, and leave the fuel in droplet form in the oxidizing 
atmosphere. 

Grouping the fuels on the basis of heat of vaporization seems 
about the only way to handle them. The other properties are too widely 
divergent. Both of the hydrocarbons, RP-1 and LCH^ can be placed in 
some sort of intermediate category based on heat of vaporization. The 
heat of vaporization of ethanol is very much larger than those of the 
hydrocarbons. Therefore, if the heat of vaporization is controlling the 
direction of spin instability, then the two hydrocarbons, having roughly 
similar heats of vaporization, would be expected to behave in a similar 
manner, which has been shown. Ethanol, having a very different heat of 
vaporization, behaves differently , as expected. These differences in 
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heat of vaporization might affect the sign of the velocity index, which 
determines the direction of spin of the wave. (The wave spins in the 
positive direction.) 

A further attempt at analysis may be made by considering 
some of the chemical properties of the fuels, specifically, their com- 
bustion traits. Not much data are available at the present time, but 
more should be forthcoming as soon as the kinetic flow reactor studies 
get underway. 

Again it is found that the hydrocarbons may be classed to- 
gether, apart from alcohol, in oxidation characteristics. The more 
fundamental data of bond strengths is first considered. All values of 
bond strength are approximate. 

k CciX 

SUBSTANCE MOLECULE BOND STRENGTH, mole ' 


ethanol 

H H 

1 1 

H-C-C-O-H 

C - 0 H 

70 


1 1 

H H 

0 - H 

95 


H 



methanol 

1 

H - C - 0 - H 

0 - H 

100 


I 

H 


H 

1 


methane 

H - C - H 

1 

C - H 

101 


1 

H 




H 

H 



ethane 

i i 

H - C - C - H 

i t 

C - H 

98 


H 

i 

H 

c - c 
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Data taken from Handbook of Chemistry and Physics , and Pauling The Nature 
of the Chemical Bond. 
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The table on the preceding page of bond strengths will at 
least give a relative idea of how the molecule will begin to react in 
a high- temperature environment. The weaker bonds will be broken first, 
and the molecule will probably react by that route rather than by another 
one requiring the breaking of a stronger bond. Also, an idea of the 
relative ease of reaction (oxidation in this case) of different molecules 
may be obtained by comparing the strengths of their bonds, and in parti- 
cular, their weakest bonds. 

Consider ethanol. The strength of its bonds are shown below 
for easy comparison: 

H H 

00] ool 
0\| cv| 

h 98 c §3 c ZP 0 ^H 

ool 001 
o>| <^| 

H H 

Ethanol will react most readily by a route which involves breaking the 
C — 0 bond and detaching the hydroxyl group. Before it loses one of 
its hydrogens, it will cleave the C — C bond and break the molecule up 
further, forming additional radicals. 

Now consider ethane. In the absence of more data for more 
complex hydrocarbons, ethane can at least give an idea of the relative 
strengths of bonds in higher hydrocarbons like RP-1. (Note that RP-1 
has no aromatic components, so it is proper to use only aliphatic hydro- 
carbons as models.) Before ethane loses a hydrogen, it will cleave the 
C — C bond. So we can expect higher hydrocarbons to react through methyl 
radicals and other small units, rather than through larger radicals. 

Methane has only one type of bond to break, the C — H bond and 
must react by some route involving methyl radicals. 

Now consider the relative rates of reaction of these three 
molecules. Clearly, the ethanol will react more quickly than the hydro- 
carbons because the OH group detaches very easily. Much more energy is 
required in the case of a hydrocarbon to cleave the C — C bond. And 
methane requires very much more than that to react. Thus it would be 
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expected that ethanol would begin to react much sooner than would 
hydrocarbons , and that methane would wait a long time before beginning 
reaction while it absorbed the energy necessary to dissociate. 

Preliminary studies in the flow reactor confirm these sus- 

21 

picions • Tests with ethanol and isopentane established that ethanol 
began reaction much sooner and continued at a fairly steady rate. After 
a significant delay, called the induction period, the isopentane then 
began to react. Its reaction occurred at a much faster rate than that 
of the alcohol once it began, though. Methane has been found to have 
even a longer induction period than isopentane. The reaction may be 
schematically pictured as: 



OXIDATION OF THREE FUELS 

— Reaction of C^H^OH - 0^ 

— Reaction of GH 3 CH 2 CH(CH 3 )CH 3 - 0 2 

Reaction of CH, - 0~ 

— 4 2 
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From these reaction rate data, together with the information 
on vaporization rate, some tentative, and probably very tenuous, con- 
clusions can be drawn about the combustion behavior of hydrocarbon and 
alcohol fuels in a rocket motor. 

It seems probable that ignition of alcohol occurs immediately 
upon injection, very close to the injector face, with a slow combustion 
rate because of the slow rate of vaporization and the slow kinetics. 
Hydrocarbon fuels, like RP-1 and methane, however, probably must travel 
some distance into the chamber before they can ignite. And then, 
because of their rapid kinetics, once the induction period has been 
passed, and their easy vaporization, they probably burn very rapidly 
in a very short distance. As a contrast, the alcohol combustion, because 
of its slowness, is probably distributed over a longer distance. So 
the combustion of hydrocarbons is probably much more concentrated than 
that of alcohol, the alcohol combustion being much more evenly distributed} 
but the ignition of the alcohol probably occurs first. 

These tentative thoughts on the location of ignition and combustion 
of the two types of fuels, together with another bit of experimental data, 
provide a vague indication of the relationship of the combustion process 
to the spin direction. 

20 

It has been found that the spinning form of instability in 
the LOX/ALC propellant system is prevented by the presence of a baffle 
attached to the injector face. A destructable baffle was used. Stable 
operation existed until the baffle burned out, and then unstable opera- 
tion began spontaneously in the spinning form. Thus, with the baffle 
present, LOX/ALC could not generate spinning instability; without the 
baffle, it did. Contrast this behavior with that observed in the LOX/LCH^ 
tests quoted above. Here, a baffle could not prevent the spontaneous 
formation of the spinning instability. Even in those tests with a baffle 
present, the spinning mode occurred. Amplitudes were significantly less, 
and some damping occurred with the baffle, but the spin was still there, 
whereas it had not been present at all for the LOX/ALC system. Now it 
is known that a baffle must have a certain minimum length in order to be 
effective, indicating that whatever it affects exists primarily in the 
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early combustion zone. But how early is "early"? The early combustion 
zone for alcohol might, as shown above, lie closer to the injector face 
than that of a hydrocarbon fuel, particularly methane. If the coupling 
which produces the spinning instability, and which determines the direc 
tion, operates near the point of ignition, it would be overcome by a 
baffle in the case of alcohol, but not in the case of methane, for 
then the baffle would not be long enough (2- inch baffles were used). 
Clearly here is an easy point to check with a longer baffle. If this 
analysis is correct, it shows that the induction time for the reactions 
- a kinetic effect - does have a significant bearing on combustion 
instability . 

Of course, these conclusions are based on very meager data, 
and so they should not be given too much weight unless they are corrob- 
orated by further tests. But the direction they point is promising. 
More testing is underway to carry the study further. 
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V. EXTENSIONS TO THE THEORY 

A. THEORY OF NONLINEAR TRANSVERSE INSTABILITY 
1. Gas Dynamic Equations 

Theoretical research work in the area of nonlinear transverse 
instability using a droplet evaporation model is continuing. The analysis 
is applicable to a thin annular chamber in which the oscillations can be 
assumed to be purely transversal - the dependence on the radial coordinate 
can be disregarded and the relevant independent variables are t , x and 
y (the circumferential coordinate) Fig. 31. For the assumed geometry the 
combustion chamber conservation equations can be derived. The following 
fundamental assumptions are made: 

1) The gaseous material contained in the chamber is in the form 
of gases of perfect combustion. The liquid, its vapors and 
intermediate products of combustion occupy negligible volume. 
This assumption is justified by the fact that observed 
oscillation frequencies are quite close to theoretically 
predicted values which implies that throughout the chamber 
the actual sonic velocity is nearly the same as that of the 
combustion gases at their adiabatic flame temperature. This 
would not be true if considerable volumes of gas were at lower 
temperatures. 

2) The combustion gases are assumed to be thermally and calorically 
perfect, and homo- compositional. Viscous and heat conduction 
effects are neglected, except for the possibility of existence 
of shock waves and the necessary heat exchange between liquid 
propellents and the hot combustion gases. 

3) The liquid phase is well dispersed in the gaseous phase and can 
be characterised by the density ^ » mass of liquid per unit 
volume of gas . 

4) The drag of the droplets is considered negligible. Therefore 
the droplet velocity is constant in space and time. 

5) The transformation into burned gases immediately follows the 
evaporation of the propellants, so that the rate of combustion 
coincides with the rate of evaporation. 
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Then the non-dimensional conservation equations for the 
chamber gases are: 
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All quantities are normalized according to a scheme wherein 
the reference state of the gas is its state at the injector end under 
steady-state conditions. The reference velocity is the sonic velocity in 
the reference state. The reference length is the circumferential develop- 
ment of the annulus. The reference time is the ratio of the reference 
length and the reference velocity. 

The RHS of the three conservation equations represent the mass, 

momentum and energy sources corresponding to the combustion rate of the 

liquid propellants. The energy source term Q coincides with the mass 

source term only because in the non-dimensionalization the energy content 

of any mass generated is taken to be unity. #/ 
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The boundary condition on the chamber walls is q^ = 0 under 
the inviscid assumption. Also at the injector x = 0 , u = 0 . To 
complete the boundary conditions the condition at the choked nozzle throat 
should be also added. Under the assumption that the combustion is com- 
pleted at the nozzle entrance, i.e. Q - 0 within the nozzle, the nozzle 
behavior can be treated as a pure gasdynamical problem and the boundary 
condition expressed at the nozzle entrance. The difficulty of solving the 
nonlinear nozzle problem can be circumvented by considering a rather un- 
conventional nozzle geometry consisting of a large number of short orifices 
A quasi-steady nozzle condition, uniformly distributed at the nozzle end, 
may be obtained by this technique. Then the boundary conditions on the 
problem can be written as: 
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local sonic velocity. 
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We assume a spinning type of instability, travelling in the 
positive y direction. Then all the quantities must be functions of 

= ft-y , where f is the frequency of the oscillations (still unknown 
at the moment) and must be periodic in x with a period of unity 

The equations can now be written in terms of cX and x . We get: 
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Besides the X boundary conditions (1,8, 1.9) we have to add the 
condition of periodicity in c< .In the case where the spinning wave 
includes a shock, the conservation conditions through the shock have also 
be be satisfied. 

The steady- state equations may also be obtained directly from the 
conservation equations. These are found to be: 




(1.14) 


t ~ / - n> u- ( u x - u) (i.i5) 

r 

h ~ / - i— u (1.16) 

— — 2 

r 

where the bar indicates 


the steady- state quantities. 

Rather than seeking the general solution of the above equations we 
restrict ourselves to the specific objective of looking for periodic solutions 
to the equations. The method followed is that of expansion in powers of a 
small parameter £ suitably related to the perturbation amplitude of the 
variables . 

2, The Droplet Evaporation Model 

The source strength Q is evaluated by assuming a model of droplet 

evaporation and combustion. An earlier analysis of nonlinear transverse 

22 

instability by Crocco was based on the n - X model (or sensitive time lag 
model) and gave promising results. However, as Crocco himself points out 
the n - T model is unsatisfactory because it fails to fulfill the ultimate 
aim of a combustion model, i.e. , the establishment of a link between the 
combustion rate and the physico-chemical processes involved in the combustion 
of the injected propellants. The droplet evaporation model is more satis- 
factory from this viewpoint. 

2 : 

Using the Priem-Heidmann model of quasi-steady droplet evaporation 
the evaporation rate is determined by the relation 


9 
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where 


K 

droplets 

Re 

Pr 


=5 dimensional quantities 

= constant related to properties of gases and of the 


Reynold’s number = 
Prandtl number 



For simplicity, at large Reynold number values the conduction 
term can be neglected in comparision to the convection term. 

Defining the relative droplet volume V = /+\ , = droplet 

radius at injection, and non- dimens iona Using the other quantities in the 
usual fashion we get ^ 

- fK ■+ “xK - ■- c (Vv) < 2 -!> 
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Further the non-dimensional evaporation rate can be shown to be 

„ D,V 

Q = -/?• - where P u, Jjl the (2.2) 

non-dimensional injection flux. 

Now in steady state we will assume the combustion to be properly 
distributed along the chamber even for vanishingly small values of the 
expansion parameter. As V = 0 this implies that = 0 also. Also 

- o(0 , p = 0^) and = Then from (2.1), for steady- state, 

u i - c ( u t) /i ° r c m °( u i) /z ■ 

We further try to determine the proper ordering of the terms. 

Consider u , u^ and Q = 0(p.) . Now v 1 , ^>‘ , p 1 , must be 

of the same order, since we want the lowest order perturbations to satisfy 

the wave equations. Let v* , p f , and p 1 = 0( €r ). We have to determine 

a relation between £ and Lu . 

22 

Crocco has considered the case when the transverse velocity 
perturbation is much smaller than . and u . The present analysis con- 
siders the case where the transverse velocity perturbation is larger than 
u^ — u = 0(Jk ) , The other case, which is not covered by this report, and 
on which work is now in progress, is when v and u^ — u are of the same 
order. 
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Ordering of the Terms for the Case v '» * 
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Therefore we have the following expansions: 
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In order to avoid non-uniformly valid solutions in the per- 
turbation approach for nonlinear problems the technique of coordinate 
stretching has been found to be very useful. We introduce here the 
stretching of both the variables o( and x in order to show the 
general method for problems of this type. (Later we find that for the case 
it is not necessary to stretch the x variable to obtain 
the solution for the waveform,) 
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The stretching expansions are: 

8 — * + ee, (e, t) -*• (e^t) -#• (2.5) 

? - •* *r t ( 9f t) t faij'*"— 

where the functions 0 f 0, f ^ etc are to 

be determined. All dependent variables are now functions of 0,1 and 
they must be periodic in 0 with period 1 . If we have a travelling 
transverse shock, the stretching of x enables us to have the shock 
correspond to 0*0, and because of periodicity the shock again appears 
at 6 * / and further at all integral values of © . 
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Introducing the expansions 2.4, 2.5 and 2.6 in the equations 
1*10 “ 1*13 and separating terms of various orders of £ we get the 
following set of equations. The equations of 0(6) are 
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Consider the £ order equations. The unsteady equation 
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as A sr O , C must be a constant 

r 'l 


Also 


ti/ Y - <P(e) + c 
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Therefore the £ order equations have the solution 
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Using the £ order equations to simplify the higher order equations 
we get. 
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Using the steady- state results, the equations for the unsteady 
c omp onent s bee ome s , 


K 


2 6 


V. 


/ 

V~ . 

26 

Y 

_ t> 


26 


26 


= Q + [( r +')’t' + rc - f,]t 6 
= - <f> 9 e l? 

= " [ -f, + c ' 4 J ’’e 


r' + c r ? = O 

26 > '? >e 


-5' 


( 2 . 9 ) 


( 2 . 10 ) 


( 2 . 11 ) 



103 


Our objective is to determine an equation for the wave shape 
<p(e) . To complete the solution we need to determine an expression 
for , the source term. 
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From the steady- state equation 1.14 
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Using this expression in 2.11, we have 
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Integrate the equation from J «■ O to / to find the average 
over the length of the combustion chamber 
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which is the equation for the wave shape. The 
integral part of the second term is of course a constant depending only 
on the steady- state solutions. Denoting it by ’iJ we have the equation 
for the wave shape ^ 

(r-n)i (0 - 1^1 Zy ~ 0 

For no value of l) can this equation provide a periodic 
solution which also satisfies the shock condition j (b clB — O . 
Integrating the above equation with respect to Q from O to I we 
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which has the 


Therefore any wave with the square root driving term (i.e. 
case for \T ' tends to decay to zero. Of course as soon as 

the perturbation amplitude decreases somewhat the square root driving 
term is no longer valid and a linear driving term has to be used. 

(The details of this case are still being worked out . For iT W U^- 
it is found that the driving term can be approximated by a linear 
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function). Thus the 3 ranges of V 1 and Mj" ** considered do present us 
with a reasonable picture of the phenomena. For \T ' <<s - U the 

dependence of the driving term (the combustion rate) on V ' is quadratic. 

22 

As shown by Crocco instability can exist (i.e., periodic solutions can 
be found). These periodic solutions constitute a trigerring limit for 
instability - any perturbation ultimately leading to an indefinite ampli- 
fication. Of course this is quite unrealistic and we should expect to 
reach a limit cycle. The reason of course is that as the amplitude 
increases the quadratic driving term is no longer valid - for large 
amplitudes we have seen that the driving term only increases as the 
square root of the V"' and leads to decaying solutions. Thus the limit 
cycle mentioned seems to exist for amplitudes lower than where the square 
root dependence holds. Thus the third range ( V~' & U^-U ) where the 
driving term is linear in v is expected to provide a stable limit cycle 
solution. As mentioned before work on this case is progressing rapidly 
and is expected to be completed in the near future. 
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B. NEW EXPANSION TECHNIQUE FOR THE CASE v'<< u - u^ 

The purpose of this article is to discuss the nonlinear insta- 
bility in an annular combustion chamber for the case when the transverse 
velocity is much smaller than the difference of liquid drop and gas 
velocities using simplified vaporization model. Two objectives will be 
sought: 1) The final oscillatory conditions; 2) the triggering limit, 

above which the amplitude of oscillations becomes amplified and below 
which it decays . 

Let us start with the gas- dynamic equations derived in Section 
A of this chapter and reproduced below for convenience. 
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-1 — 6 “ 

with p __ >|p & Q * 

this set of equations describes the flow phenomena in an annular rocket 
combustor. Discussion of the assumptions for these equations will not 
be given here. (See Part A). With the short nozzle approximation which 
gives the boundary condition at the end of the combustor, our main objective 
is to seek a periodic solution for spinning waves. At the nozzle end, the 
boundary condition is 
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The steady- state solution can be obtained simply by integration: 
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The injection rate is expressed by 



The source term Q is being evaluated by assuming a vaporization 

model, i.e., the combustion process immediately follows the propellant 

25 

vaporization. A simplified model has been suggested (also see Part A) 

q = c ( vp % r <«> 

and the equation for the conservation of liquid drop mass is: 


\Jc( u* V*, •= — Q 
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where V is the drop volume. 
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The method of solving this problem is that of expansion in 
powers series in a small parameter £ , which is related to the 

perturbation amplitude of the physical quantities. After introducing such 
expansions in Equations (1) to (4) and (11), we obtain equations which 
correspond to different powers of 6 and can be solved separately and 
successively. For small amplitude and vanishing gas source strength the 
steady- state conditions are uniform and the perturbations should satisfy 
the acoustic equations of which the general solution is known in some 
particular cases. In Refs. 22, 24 and 25 we see that choice of the expansion 
parameter £ is very critical and can lead to different results according 
to the different orders of the steady and perturbation quantities. It is 
the purpose of this discussion to show a proper way of selecting such a 
parameter which is based upon some essential requirements. We should not 
be mistaken that these requirements are the only ones, we could expect 
some modifications for other cases. However, the treatment below does 
have its own generality. 

Let us assume that the unsteady parts of pressure, gas density, 
entropy, transverse velocity component and drop volume are of the same 
order, namely 


i>' = o o') 
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P' = o (V) 
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<5-'= 0 (V) 
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v'= OCV) 
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the relation between *\J and has not been established yet, but will 
be determined later. 

From the expression of the combustion source (10) where 
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and with the assumptions 
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the first relation of (17) indicates the fact that we shall consider the 
transverse mode only and (17) and (18) are the small perturbation assump- 
tions . 
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with the use of (20), From (13) and (24), we obtain 
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substituting in (13), (14), (15), (21) and (23) we have 

p ' '= oc%') 
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If we assume that the derivatives of all perturbation quantities 
with respect to ©( and ZC are °f the same order as the quantities them- 
selves, i.e. , 


-fu' = O(-ip') , t-4 = o(t') 

W = OCt/) 

since ^=OC-i) , U = Ogu) 


after separating the steady- state part from Equation (1), the left-handed 
side of the unsteady equation is of the order 

LHS = 0 

and the right-handed side is 

RHS = 0(Q’) 
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If we expect that the first order equation of the unsteady 
quantity be the acoustic equation (it should be realized here that the 
first order does not mean the term of power £ , rather, it is the first 

leading term of the unsteady quantity), clearly Q' should not appear 
at this order, the necessary condition would be that Q' is of higher 
order than v' . 

Since 6 is the actual expansion parameter, if we choose 

q' — OC'ir'e) (33) 


the above condition can be satisfied. Of course, the cases with 
Q / = 0(l/e > ')for any riZ l would serve the same purpose but obviously 
will give different results. We shall not discuss other possibilities 
in this article. 

The above argument for getting an acoustic type equation can 
also be applied to the drop volume equation. With the assumption 

Vi = 0 cv') , v* = o(v') 

we obtain the unsteady- state equation from (11).. The left-handed side 
is of order 

LHS = OC V') = OC-vV) 

and the right-handed side 

RHS = O(^Q') 

In order to obtain an acoustic type equation, we choose 

q / =~ O C'u'e ) 

From (33) and (34) , should be the same order as jj, 

U.J, = Oiji) 

this is what we observe in many practical rocket engines. Some of the 
22,25 

former analyses require that (jl^ and JLA be of different order. 
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<2 = 0 (v'e) 


it is seen that 


0(j^) 


For convenience, we take . The final result indicates 

that the expansion parameter should be the injection rate which is 
approximately the Mach number at the nozzle entrance. 

Since we have obtained a proper expansion parameter, we can 
proceed to expand every variable into a power series of JX and substitute 
into the equations* After separation of terms of different powers in , 
the equations can be solved separately and successively. In many nonlinear 
problems, it is found that method of coordination stretching is very useful. 
Besides expanding the dependent variables, we also expand -p ^ , 31 in 

power series of JX and a set of independent variables and Q is 
introduced to replace 1C and o{ . The following expansions have been 


taken 


C ■=/*-' 


l -‘rJJ, -ip, + JLX-pz 


(° — ' + y u (°i ' + pz + 




JO? -4-yU 1 

JJ< -4yU 


&-L 
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<1 

II 

53 

4yM Vi -tyU 2 V/4 -VyU 3, 14 

\7~ 

yV* Vi 4- yU* v// 

u = 

^aH.| CX-i 4y(A (-X j +-yA 

u'= 

> 

V = 

O 

u'~ 

ym 2 *uV 4y/U* Vi 

Ut = 

^xCLe, 

Q = 

^/X. <5, -+yM <j?i 4-^M? 

Q* — 

yV 3 6? a 


^ f . +/ I- - I s 




with < g, — U| - u«, 

U*. 


( 48 ) 
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The coordinate stretching 

■f — -X- ■% 3 

e = <=( + jx ©i +m 2 ©*. ®j 

-f - \ + /A+, +}*.*■$■,. -V/A?-- f s 

0 are t ^ le independent variables. 

From (1), (2) and (11), we obtain the first order equations 


(49) 

(50) 


vanish, 


U_i-| — 

(51) 

II 

0 

(52) 

Vo t C t.o -V- u«, ) = -u«, (5, 

(53) 


Since ^ C'f— O') ~ O j from (52) we conclude that *^> f should 


i.e. 


^lp. = o 


(54) 


the steady part of (53) is 


Vo t =-Q, 

the expression of Q x is 

Q, == — c Vo f, )^= C V. 

where — ^ ^ 

* ^ 

The differential equations (51) and (55), combined with (56) 
and also with the boundary conditions that ]/ 0 (_o) — d Cl, (0^=-O« yield 
the following results: 


(55) 


( 56 ) 
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For 


° < U, ^ U«, ^ o ^ f $ f 0 

Vo ■= ( C vJL, C-6r| ) — ^VcvJL C-fe-f ) 

U ‘ ~ l~ (0«A N C4t') “aJ^ MJL CM)t 


(57) 

(58) 


is determined by the condition that at ^ = o (jj^ 

For $ a, ± , 1= =s t $ 41 

V„ = ‘v 1+ U-U^-f.'fl] 




u» ■== 


^Cv+-u«, ) + o-u<, )Oo\;2-fe(-f-f„)')j 


where is the end of the combustion zone, namely, \/ e> (Q') —O , and 

u,m=\ • 

4 -- -JL- -+ 

^ 2-^ So 

The unsteady- state equation of (53) is 

^o< *^if© == O 


since \/<»-« is not identically zero, so it is necessary that 


t 


IO 


O 


(59) 


(60) 


(61) 


(62) 


( 63 ) 


that means ^ ( is function of ^ only or ^ ( ■=. . 

Since there are no unsteady terms in > we conclude 

that -f, must be zero because this is the frequency correction to that 
order 
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Now we go to the second order equations. They are; 


l -T7 - 5 - 1 




± 


* +G >' u ' =~Q. 


t, 


2© 


■u^ = 


v - ,2e = o 


<SVe -+- U, <S"j.g — - (3, O', 

C \4e> — Q, 5z e ~] *-<■«, ( — Q ^ 

r w - G*. 


— U* 


* " 2 . 


Vc 


& 


J 


(64) 


(65) 


( 66 ) 


(67) 


(68) 


The steady- state equation of (67) reads 

?v 


©i* H- (h 

l s TX, u » 


the solution is 


o 


f ^ 


- A 


u., 


but at ~%=o > — O this implies that f\ should be zero, or 

(?, •= o 

The unsteady- state equation is 

=0 


(69) 


(70) 


( 71 ) 
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the solution should be 

6^ z“ : b '=: Co^\s'[~ 

if there is any dependence on ^ at all, it can be included in „ 


(72) 


The steady part of equation (68) 


V, 


it 


Qy 


v, 


Q\ % % = - [■ 

* v Vo 




and (64) can be solved simultaneously. Rewriting these two equations, 
we have __ 

= Q, -S, 


v/, 


is 


Vo 


u, 


a, 


V \ / 


Vo 


-V- 


< &» 


■] = 


(73) 


(74) 


(75) 


Adding these two equations and integrating, we get the result 
after applying the boundary condition that V,(o) = lU(o)= 0 . 


V\ -+ -U. x 


o 


(76) 


so (74) becomes 


(}L\ _ (• 

V v. J-f 2- L 


v. + 


1(A) =-^-t 

3 \/ 0 \/o 


(77) 


The general solution is 

V, _ 


Vo 


-4“ ~C-^c^p 


a "t- 1 £ ,J I- 


■ {~r [- £’-f 1 s. I 1 * 5 
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The boundary condition requires that A - — O * 

The choice of should be such_that the power series 

ve 

» the solution would be 

V, = O 


expansion of V converges everywhere, i.e., should be finite 

we see that if we take 


Hence , 


(78) 


also 


u 


o 


(79) 


(65) yields a steady- state equation: 

Jh =:-Q, U U, - U< L, ) 


(80) 


which gives the result 

u, C U, - Ue, ) 

0 

with the boundary condition -felo) — O . 

The unsteady-state term of (65) is 

rU --=o 

From (66) and (82), we get 

A A ~ 4 

Jr 

lr/ — 4> (© y-v- c (?) 

where <£> is function of © only and C is function of only and 
they will be determined by the higher order equations. 

We now consider the third order equations: 

From (4) , we obtain 

G*/© ■*" ^ » ft $ ^ 1 ~ ^ 

- Uu*,u, - a* )3 


(81) 


(82) 


(83) 

(84) 


( 85 ) 
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the steady- state relation is 


Sl s + & = Q, (^K*R, -4U<,) 

With the boundary condition S^(o )—0 > the solution is 


= 


is- -4 , iT 


cu, -2U <( ) U, 


The unsteady- state equation is 


^4 - <3-^ H 


-K - 


<2, t (MO 4 - 


From (1), we have 




jSl / 

* "Vo 


- < 3 , 1 
+ Q,[ 


SCfr-*-0 


s-. 




1-2 lr 


i-t (X.) 4-±(l±C 

J T7 ' 4. ^ rr 


after substitution of the results from the first and second order solutions. 
The steady- state equation is 

u _SL(iL + -iL} 

* _ w > 
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and the unsteady- state equation is 



We obtain the following equation from (2) 


('Hi + -He <9, t ) + Q, == O 

The steady- state part is 

=° 

which gives 

+1=0 


( 91 ) 


(92) 


(93) 


(94) 


because of the boundary condition. 

The unsteady- state equation is 


" + ©l t + Q* ?». = 


o 


From (11), we have 


[ ( We — Q\ ~Z-ie ) \Af ^ 

+ ««,( v, s +4 + 


/ i 


U, -V- 


_ir /i* 


U, 


'M 


(95) 


(96) 




■W 
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The steady- state equation: 

v + <3i ( 4l- +• ) 

+ 2. V Vo ^ ; 


(97) 


~ ~ ^ U«, M, - %. f ^ 


and the unsteady- state 


( V *© — Qi ~? 5e ') 

= .-^,5, ^-4 c a-4) + x(^) + i ( ^cj 

We add (90) and (97) together and use the boundary condition 
that V z ( 0)-0 and ( o) = O , we obtain 


? (98) 


= (-^|L u, - Js-Cfe, ) u, x <"> 

Substituting this result into (96) and (97), we have the 
following two equations: 


( £ 

) 

<3, 

(4r+4r 

V- 


1 




Vo 

% 

/ V 

¥. ' 


— 

h 


4- -Hj 
W 

)(- 

M _ 

“ “t" ) (100) 



-e - 

^ Q, 

i ~ 

"fj 3 

1 


(£ 

Vo 

\ ~ 

2. 

- ( v. 

+ i 

')( 

1 

u 

) 


— — < 

3, 

TS 


u, 

8-. 

)(- 

“2. 

CU, U, ^ S* J (ioi) 
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The solutions of these two equations present some difficulties 
here because of the singularities at — O and |/ 0 =rO • By taking 

(no stretching) 4 straight foreward integration yields the 

2 5 v T? CL 

results that 11 and r - tend to infinity at the singular points 

and this is not acceptable because the requirement for the validity of 
series expansions, and —M* should be finite everywhere. Also the 

singularities are the^causes ^'difficulties to solve the unsteady 
equations. Although no final results can be concluded at the moment, 
we have been able to show the correct approach to this problem. It is 
believed that the technique or ordering in this article can be extended 
to other similar situations. Work has been progressing and will be 
completed soon. 
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